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THE FEASIBILITY OF A GLOBAL AIRLIFT* 


This paper is a reprint of a study prepared by R. S. 
Weinbergt and J. W. Higgins of the Naval Warfare Analysis 
Group. This group is maintained by the Massachusetts Insti- 
tute of Technology for the Long Range Objectives Group in the 
Office of the Chief of Naval Operations. 





Note by Managing Editor: The impossibility of a global airlift can be 
established by very simple means intuitively. However, the methods 
used by the authors in arriving at this same conclusion, even by a 
limited analysis, may be of interest to the reader. This paper gives 
quantitative proof of the magnitude of effort required for importing by 
air in lieu of using sea lift. 


H. P. J. 











CONCLUSIONS 


e This analysis determines the number of very large transport aircraft and amounts of 
fuel needed to airlift ali U. S. imports and exports for their level in 1953 (55.5 million tons); to 
airlift only vital strategic materials (9.6 million tons); and to operate an airlift for both these 
levels of commerce where the aircraft fuel needed to extend cargo aircraft range is transported 
by surface tankers. Finally, the number of dry cargo ships needed to provide equivalent sea 
transport is determined. The resources needed for each of these three cargo operations are 
summarized in table A, in terms of aircraft, aircraft fuel, tankers, and dry cargo ships. 

e The demand for petroleum resources and numbers of aircraft for a full global airlift 
is so great as to preclude such an operation. 

e The demand for petroleum resources and numbers of aircraft for an airlift limited to 
carrying vital strategic materials would make such an operation very costly, although within 
U. S. capabilities. 

e Sea transport of airlift fuel would require a shipping protection effort almost as great 
as that required to protect dry cargo sea transport by modern cargo ships of the Mariner 
class. 


e@ The trend toward greater dependence on overseas sources of heavy materials will 
outweigh foreseeable technological improvements in aircraft design and production, which 
might otherwise permit limited aircraft competition with surface ships. Large transport 





*Manuscript received January 8, 1958. 
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FEASIBILITY OF A GLOBAL AIRLIFT 


TABLE A 
Resources Required to Operate Each of the Three Overseas Transportation Systems 
(A Conservative Estimate Based on 1000-Miles Between Airfields) 





Mixed 
Air-Surface Lift* 


Complete 
Airlift 


Complete 
Surface Lift 





Resources Required 


Strategic 
Materials 
Airlift 


Full Lift 


Strategic 
Materials 
Airlift 


Full Lift 


Strategic 
Materials 


Airlift 





Super cargo/tanker 
aircraft 





Aircraft fuel (millions 
of long tons) 





Notional surface tankers 





Notional dry cargo ships 
or 
Mariner class cargo ships 





Naval fuel (millions of 
long tons) . 5.8 99 
to to 
9.3 1.6 





























*Aircraft and aircraft fuel values also apply to situations where fuel is available at necessary 
locations without being shipped from U. S. 


aircraft can serve as very useful supplements to sea transport but cannot provide a practical 
substitute for any substantial part of present sea transportation needs. 


I. INTRODUCTION 


A frequent proposal for continuing our foreign commerce if our control of the seas is 
strongly contested is to carry on some or all of our overseas trade by air. The feasibility of 
such a global airlift is alleged to rest in the recent development of very large aircraft which 
might even be capable of replacing most merchant-type shipping. 

This study makes only a limited economic evaluation of the feasibility of a global air- 
lift. There are, of course, other considerations which could affect the practicability of such a 
large airlift. For instance, an enemy interested in destroying our commerce would undoubtedly 
divert some of his anti-shipping investment to systems for attacking the airlift. As a result, 
large investments on our part might be required to defend the transport air bases and even the 
aircraft routes. A large-scale airlift will also require many skilled air crews, air bases, 
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maintenance facilities, and other support organizations and equipments. None of these require- 
ments are considered here. 

Because aircraft fuel and aircraft themselves represent the resources which are likely 
to be in shortest supply in wartime, only the demands for cargo aircraft, tanker aircraft, and 
their fuel are analyzed here. This is done for two types of airlift: a complete global airlift of 
all foreign commerce, and a partial airlift of only the most vital minerals. In each case the 
requirements are derived for the complete airlift of the cargo, and of the fuel which must be 
flown out for the return of the cargo aircraft. In addition, the number of "notional" surface 
tankers which would be needed to substitute for the aircraft tankers is determined for both the 
full-scale global airlift and the limited airlift. Finally, the requirements for cargo vessels and 
ship's fuel are determined for both the full-scale overseas lift and the limited lift of only vital 
mineral imports. 


Il. METHOD OF ANALYSIS 


Consumption of aircraft fuel and required numbers of aircraft can be determined as a 
function of the total distance to be flown and the distances between adjacent fields on the airlift 
routes, given estimated operating characteristics of the cargo and tanker aircraft to be used. 
The hypothetical turbo-prop aircraft considered here is of somewhat higher performance than 
any yet designed; it has an estimated capability of transporting 120 long tons of cargo, minus 
fuel used in flight, at an average speed of 400 knots and with a fuel consumption rate averaging 
-015 long tons per mile. It is assumed that these aircraft can be utilized on the airlift for 200 
hours per month. 

Table A-1, appendix A, shows fuel requirements to lift one ton of cargo for distances 
ranging from 500 to 10,000 miles, under various assumptions as to length of each step of the 
airlift and availability of fuel enroute. From this table it is easy to compute specific require- 
ments for aircraft and fuel for given tonnages to be airlifted. 


I. FULL SCALE AIRLIFT 


Table I shows our total 1953 imports, as a function of the distance they were carried, 
while figure 1 shows the principal routes over wnich these imports might be carried by aircraft. 

The tonnages in table I were based on import and export data for 1953 prepared by the 
Department of Commerce. From these total imports and exports, all trade with Canada and 
Mexico was subtracted, since this could be moved overland or by the Great Lakes. All 
petroleum imports and exports were also removed, because these are not essential to survival 
so long as our domestic resources are adequate. Imports from and exports to a given overseas 
point were then matched where possible, so the aircraft would be loaded in both directions. 
Where imports and exports were not equal, the lift requirement for the larger of the two was 
the figure used in table I. 

With the aircraft and fuel requirements per ton, and the total tonnage distributions in 
table I, the total amount of fuel and the total number of aircraft can be computed. These are 
shown in table I. 

Table II clearly shows that the cost of a complete global airlift would be prohibitive. 
Actual fuel and aircraft requirements would lie somewhere between the 1000-mile and 2000- 
mile values, since trans-Atlantic and trans-Pacific airlift routes all must involve individual 
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TABLE I 
Distribution of the United States Imports and/or Exports According to 
the Estimated Airlift Distance Between Their Points of Origin and/or 
Destination and CONUS 1953 





Airlift Distance Payload Airlift Distance Payload | 
(miles) (10° long tons) (miles) (10° long tons) 





500 6225 5500 5750 
1000 3775 6000 2950 
1500 1850 6500 925 
2000 1075 7000 1250 
2500 1075 7500 2050 
3000 4400 8000 2150 
3500 7875 8500 1150 
4000 5425 9000 275 
4500 2600 9500 50 
5000 4655 10,000 _ 
Total 55,500,000 Long Tons 




















Total Cargo Movement = 221,325 million ton miles in the direction of heavier 
traffic (individually for each route) 


TABLE 
Fuel and Aircraft Requirements to Airlift a Flow of 
Foreign Commerce at a 1953 Level 
(Aircraft requirements make no allowance for aircraft attrition) 





Airlift with Airlift with 
2000-mile steps | 1000-mile steps 





Fuel Requirements: millions of long tons of 
JP-4 type fuel, per year 


Aircraft Requirements: number of hypothetical 
aircraft assigned to airlift 16,010 10,950 

















steps in the 1300 to 2000 mile range. Even if airfields (and therefore suitable islands) were 
spaced at 1000-mile intervals, a full airlift of normal foreign trade and supporting fuel would 
require fuel equivalent to 46 percent of our 1956 petroleum production, or about 43 percent of 
our 1956 petroleum capacity. An airlift of this size would exhaust our 1956 known U. S. oil 
reserves in about 30 years. The 11,000 very large aircraft would represent 1850 million 
pounds of airframe, or a production program nearly 75 percent of our total World War II 
airframe production. 
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IV. STRATEGIC IMPORTS COMPREHENSIVE AIRLIFT 


The results of section III show quite plainly that an airlift could not possibly carry all 
of our foreign commerce. This section will, therefore, consider the feasibility of airlifts 
which are restricted to only vital mineral imports. Table III lists 30 strategic minerals, with 
the quantities imported in 1953. If imports from Canada and Mexico are deleted from this list, 
and if only extremely limited amounts of special iron ore are air-lifted, the annual require- 
ments in table IV are generated. 

The annual requirements tabulated in table IV are quite conservative. For instance, 
no organic materials are included although some of these may assume a vital position. Fur- 
thermore, no substantial quantities of iron ore have been included in table IV (although listed 
in table III), because Canadian sources, present domestic mines, and the growing exploitation 
of domestic taconite may satisfy demands. If this were not true, our airlift needs would be 
substantially multiplied. On the other hand, our most vital organic requirements could 
probably be met by a 10 percent increase in the total tonnage considered in this limited airlift. 

Following the same procedures as those in section III, the requirements for fuel and 
total number of aircraft can be derived, using the fuel consumption rates per ton in table A-1. 
These results are given in table V under the assumption that exports to the supplying countries 
are airlifted and are equal to imports in weight, and alternatively that no U. S. exports are 


TABLE Ii 
Annual Imports of Selected Strategic Materials - 1953 
(long tons) 





Complete Partial 





Strontium minerals 8800 Cadmium 700 
Columbium-Tantalum 500 Cobalt 9800 
Beryllium 7400 Copper 600,000 
Chromite 1,964,300 Iron Ore 10,700,000 
Manganese 2,946,700 Lead 270,000 
Nickel 116,400 Mercury 2900 
Platinum 300 Selenium * 
Tin 110,000 Tungsten 27,700 
Asbestos 628 ,800 Zinc 610,000 
Corundum 5000 Fluorspar 323,000 
Industrial Diamonds * Iodine 400 
Graphite (strategic) 5100 Aluminum Ore 4,480,000 
Jewel Bearings 200 Antimony 17,800 
Mica (strategic) 7000 Bismuth 300 
Quartz (radiograde) * Titanium (Rutile) 15,000 

Total 5,800,500 17,057,600 
5,800,500 
Total 22,858,100 














*Relatively trivial (Under 10 long tons) 
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TABLE IV 
Distribution of United States Strategic Imports by 
Distance from Origin to U.S. - 1953 





Airlift Distance Payload Airlift Distance Payload 
(miles) (thousands of long tons) (miles) (thousands of long tons) 





500 ---- 5500 441 
1000 2500 6000 768 
1500 1515 6500 436 
2000 37 7000 593 
2500 156 7500 1021 
3000 392 8000 519 
3500 358 8500 10 
4000 316 ---- ---- 
4500 309 ---- ---- 
5000 222 ---- ---- 
Total 























Total Cargo Movement = 37,343 million ton miles 


airlifted, so that the full capacity of the airlift is available for transporting its fuel from the 
U. S. to airfields along the route. The computations are based on availability of aircraft fuel 
without airlift at all U. S. bases and at Dhahran, the 6500-mile point on the 7500 to 8000-mile 
airlift from Pakistan and India. This route and the Malaya and Indonesia route (via Hawaii) 
are the only ones in the over 7500-mile bracket. The India strategic import tonnage is 
dominant, since India is America's principal supplier of manganese. Availability of fuel in 
Venezuela would reduce fuel consumption by about 10 percent on the South Atlantic airlift but 
would not affect Caribbean imports (Jamaican bauxite, principally) since southbound aircraft 
would have more than enough capacity for their own return fuel. Indonesian oil is assumed not 
available, and since imports from Indonesia and Malaya are small in tonnage anyway, this does 
not affect the conclusions as to feasibility of a comprehensive strategic materials airlift. 
Table V shows that the strategic minerals airlift will still require more than 12 million 
long tons of fuel per year, while over twice this amount will be required if an equal weight of 
exports must be airlifted. The minimum figure of 12 million tons is not quite 4 percent of 
1956 petroleum production, while the minimum aircraft requirement of 822 is over 10 percent 
of the U. S. aircraft industry's peak annual World War II airframe weight production. The 
higher (export and import) figures correspond to 8-10 percent of our 1956 petroleum produc- 
tion, and to 20-30 percent of our peak annual World War II airframe production. This reduced 
airlift is well within the capability of American industry but is still a large and costly operation. 


V. MANGANESE AIRLIFT 


Before considering such alternatives as sea transport of aircraft fuel, requirements 
for airlifting manganese are presented in detail as a special case. Manganese is a vital 
strategic material and is imported in substantial amounts from distant sources. The analysis 
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TABLE V 
Fuel and Aircraft to Airlift Strategic Imports at 1953 Levels 




















Exports Airlifted Airlift Limited 
to Balance Imports to Imports 
2000 mile | 1000 mile 2000 mile | 1000 mile 
steps steps steps steps 
Fuel Requirements: millions of long 
tons of JP-4 type fuel, per year 33.0 23.5 13.7 11.8 
Aircraft Requirements: number of 
hypothetical-type aircraft assigned to 
airlift 2291 1630 953 822 























shows that airlift requirements are very sensitive to sources and routes assumed, and that 
development of nearby sources of supply, in this case Brazil, can reduce fuel and aircraft 
needs to a level which could be met without strain on our resources. The minimum require- 
ment, however, still includes a ton of fuel consumed for every ton of manganese imported. 

During 1953, the United States imported nearly three million long tons of manganese, 
much of it from India and Africa. More recently development of manganese resources in 
Brazil has been intensified. A 1952 import breakdown from Bureau of Mines data lists 
2,485,900 long tons produced in or imported into the U. S. during that year; of this, 2,382,100 
tons were imports (1953 imports were 2,946,700 tons, as in table III) and 2,029,500 tons came 
from outside the U. S., Mexico and Cuba. 

Table VI summarizes the fuel and aircraft requirements to airlift 2,029,500 tons of 
manganese to the United States from our 1952 sources of supply, using the Dhahran route and 
Saudi Arabian oil for the Indian manganese airlift, or alternatively bringing manganese from 
India via Manila and Hawaii. A third possibility is also considered; that of meeting our airlift 
manganese needs entirely in Brazil. These computations, more detailed than those of sections 
If and IV, include the availability of Venezuelan oil at Trinidad on the South Atlantic and 
Brazil routes. The comparison clearly shows the sensitivity of airlift requirements to routes 
and distances, and also shows how airlift requirements can be reduced to a practicable level 
by developing nearer sources of supply, even if these sources cannot be linked to the U. S. by 
land transportation. 

While transoceanic airlift costs will be 15 percent to 60 percent higher if 2000-mile 
steps are used, the Brazil airlift is less sensitive (14 to 23 percent increase) and a route to 
Brazil with airports each 1000 miles or less should be relatively easy to establish. The 
Dhahran route can be operated with airfields at 1000-mile intervals except for the 2000-mile 
U. S. - Azores segment. Distances between airfields on the trans-Pacific route (the "Dhahran 
closed" alternate) can be held to 1000 to 1400 miles except for the 2080 miles from San 
Francisco to Oahu. 
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TABLE VI 
Fuel and Airlift Requirements to Airlift Manganese Imports 
at 1952 Level 





Dhahran Open Dhahran Closed From Brazil Only 





Total Tonnage 
Airlifted to U. S. | 2,029,500 Long Tons | 2,029,500 Long Tons | 2,029,500 Long Tons 


Fuel Requirements, 
1000- Mile Steps: 
If only imports 


are airlifted 4,009,300 Long Tons 5,968,600 Long Tons | 2,009,200 Long Tons 
If an equal weight 
is exported 6,675,600 Long Tons | 15,297,600 Long Tons | 2,642,400 Long Tons 


Aircraft Needed, 
1000- Mile Steps: 
If only imports 


are airlifted 28 42 14 
If an equal weight 
is exported 47 107 19 




















VI. COMBINATION SURFACE AND AIRLIFT 


So far, foreign trade has been considered as an exclusive airlift. Since the preceding 
sections have shown that a complete airlift is impracticable, and that transoceanic airlifts of 
strategic materials, though conceivable, are costly, one compromise solution might be the 
airlifting of cargo with a surface lift of the aviation fuel needed for the cargo aircraft. Since 
aviation fuel requirements are so enormous for both full and strategic imports airlift, very 
substantial economies might be effected by using surface tankers to supply this fuel. Further- 
more, tankers are among the fastest and most modern segment of the U. S. merchant fleet and 
should therefore be easier to protect from enemy action than much of our cargo fleet. Table 
VII shows the distribution of the fuel needed for cargo alone for both airlifts, indicating where 
this fuel must be made available. 

From table VII the full-scale airlift can be seen to need 63 to 72 million tons of fuel 
annually, of which 56 million tons must be moved to airfields outside the U.S. The comparable 
figure for an airlift which provides its own support is 158 to 231 million tons of fuel, depending 
on the step length—which shows that 60 percent to 68 percent of the aircraft fuel requirement 
can be eliminated by a sea lift of enough fuel for cargo aircraft. Comparable figures for the 
strategic materials airlift, with exports airlifted to balance imports, are 10.6 to 12 million 
tons of aircraft fuel, of which almost 9 million tons must be moved from the U. S. or from 
Dhahran to other areas not linked by pipeline. If the strategic imports airlift provides its own 
support, it will require 23.5 to 33 million tons of fuel, depending on the step length, so that the 
fuel sea lift saves 55 percent to 64 percent. 


506275 O-59 -2 
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TABLE VII, Part A 
Cargo Aircraft Fuel Requirements 


Distribution of Ocean-Transported Fuel to Support Full-Scale Cargo Airlift 





Distance from the U.S. | Airlift with 2000-Mile Steps | Airlift with 1000-Mile Steps 
(miles) (thousands of long tons) (thousands of long tons) 





0 (at U. S. bases) 15,921 7464 
500 417 417 
1000 540 13,414 
1500 427 124 
2000 26,320 12,250 
2500 72 72 
3000 629 10,643 
3500 1819 530 
4000 15,509 7336 
4500 174 174 
5000 666 4962 
5500 1328 385 
6000 5521 2597 
6500 62 62 
7000 179 1646 
7500 474 137 
8000 1335 642 
8500 77 TT 
9000 39 50 
9500 12 3 











Total required: 72,021 63,035 
Total long tons to be 

transported: 56,100 55,571 
Million ton-miles of fuel 

transport: 187,100 164,300 





Aircraft required for airlift: 5000 4380 








If the fuel required (table VII) is carried by T2-SE-A1 tankers which make 14.5 knots, 
carry 130,000 barrels (15,800 long tons) of aircraft fuel, require 8 days in port per round trip 
and have a repair requirement of 13 percent of cruise time, tanker requirements for this fuel 
lift will be as shown in table VIII. 

Operating the surface tankers will require 2.3 to 2.6 million long tons of naval fuel for 
the full-scale airlift, and .33 to .39 million tons for the limited airlift. 

The advantages of the mixed surface-airlift are obvious: in the full-scale lift, less than 
300 notional tankers and about 2.5 million tons of naval fuel per year can replace 6600 to 11,000 
large aircraft and from 95 to 156 million tons of aircraft fuel. In the limited airlift, 40 to 44 
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TABLE VII, Part B 


Cargo Aircraft Fuel Requirements 


Distribution of Ocean-Transported Fuel to Support Strategic Imports Airlift; 
Cargo Aircraft Loaded in Both Directions and do not Double as Tanker Aircraft. 





Distance from the U.S. 
(miles) 


Airlift with 2000-Mile Steps 
(thousands of long tons) 


Airlift with 1000-Mile Steps 
(thousands of long tons) 





0 (at U. S. bases) 

500 
1000 
1500 
2000 
2500 
3000 
3500 
4000 
4500 
5000 
5500 
6000 
6500 
7000 
7500 
8000 
8500 


Distance from 
Dhahran 
0 (at Dhahran) 
500 
1000 
1500 
1000 (by pipe- 
line) 





2565 
0 
358 
350 
3310 
10 
935 
83 
1962 
743 
32 
102 
893 
29 
85 

0 

91 

1 





1372 
0 
2271 
102 
1578 
10 
1320 
293 
891 
387 
716 
30 
448 
29 
162 
0 

39 

1 








Total fuel required: 


Total long tons to be 
transported by sea: 


Million ton-miles of 
fuel transport: 


Aircraft required for 
cargo airlift: 














TABLE VIII 
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"Notional" Tankers Needed to Surface-Lift Cargo Aircraft Fuel 


























~ 
a | @ | @ (4) (5) © mM | @  @ | 
Tanker | Trips Annual Tankers | Tankers Needed to| Tankers Needed to | 
Lift Round-| Per | Delivery Rate | Needed for} Support Full-Scale | Support Strategic | 
Distance] Trip Year Per Annual Lift Cargo Airlift: Imports Airlift: 
(miles) | Time Per | Tanker (000 | 1,000,000 | 2000-Mi | 1000-Mi | 2000-Mi | 1000-Mi | 
(days) | Tanker LT Fuel) LT Fuel Steps Steps Steps Steps 
500 11.2 32.6 515 1.94 8 8 -00 .00 
1000 14.5 25.2 398 2.51 1.4 33.7 .90 5.70 
1500 417.7 20.6 325 3.08 1.3 4 1.08 31 
2000 21.0 17.4 275 3.64 97.6 44.6 12.05 5.74 
2500 24.2 15.1 239 4.18 oa 3 .04 .04 
3000 27.5 13.3 210 4.76 3.0 50.7 2.55. 6.28 | 
3500 30.7 11.9 188 5.32 9.7 2.8 44 1.56 | 
4000 34.0 10.7 169 5.92 91.8 43.7 11.62 5.27 | 
4500 37.2 9.31 155 6.45 eS | 1.1 4.79 2.50 
5000 40.5 9.01 142 7.04 4.7 34.9 23 3.04 
5500 43.7 8.35 132 7.58 10.1 2.9 ott 26 
6000 47.0 7.77 123 8.13 44.9 21.1 7.26 3.64 
6500 50.2 Vom 115 8.70 a) Dd °29 -20 
7000 53.5 6.82 108 9.26 4a 15.2 -78 1.50 
7500 56.7 6.44 102 9.80 4.6 1.3 .00 -00 
8000 60.0 6.08 96.1 10.4 13.9 6.7 .09 .04 
8500 63.2 5.78 91.3 11.0 eB a 01 01 
9000 66.4 5.50 86.9 11.5 0 a .00 -00 
9500 69.7 5.24 82.8 12.1 0 0 -00 .00 
Total Tanker Requirement: 288 261 44% 40* 
Fuel needed for tankers, thousand long tons 2600 2300 390* 330% 























*(Includes one tanker based at Dhahran and not included in the body of columns’8 and 9.) 


(Notes on computation of table VIII: 


Col 
Col 
Col 
Col 
Cols. (6) through (9) apply col. (5) to the 


» (3) = 





. (2) = 1.13 [2 x Col. (1)/14.5 x 24] + 8.0. 
= 365 + Col. (2). 
. (4) = Col. (3) x 15,800. 

. (5) = 1,000,000 = Col. (4). 
requirements given in table VII). 
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notional tankers, using less than half a million tons of naval fuel, replace 900 to 1460 large 
aircraft and 13 to 19 million tons of aircraft fuel. Sea transport of fuel reduces total airlift 
requirements by over 50 percent. 


Vil. COMPLETE SURFACE LIFT 


As a means of final comparison, this section develops the requirements for surface 
lifting all of the cargo in both the full-scale and the limited strategic overseas trade. To lift 
this cargo a notional dry cargo ship has been used with these "average" characteristics: 
sustained speed, 10 knots; average payload, 7900 tons; port and repair time per round trip, 

39 days. Alternative assumptions use ships of the Mariner (C4-S-1a) type for the full-scale 
foreign trade lift and either Mariners or C5-S-AX1 bulk cargo carriers for strategic mate- 
rials. Basic capabilities are 10,000 tons at 20 knots for the Mariner and 22,000 tons at 16 
knots for the C5 bulk carrier. 39 days port and repair time per round trip is used in all 
cases. With these assumptions the numbers of cargo ships needed can be computed as a 
function of distance, as in table IX. 

The cargo ship requirements in table [IX can be applied to the cargo requirements of 
tables I and IV to derive the total numbers of ships needed for the surface lift. The full-scale 
lift will require about 1390 notional 10-knot dry cargo vessels and about six million tons of 
bunker oil, as compared to over 10,000 large aircraft and over 150 million tons of aircraft 
fuel for the full airlift. 846 Mariners will also meet the full-scale lift requirements, increasing 
fuel needs to about 9 million tons of bunker oil. 

The limited overseas lift would need about 238 10-knot ships, 145 Mariners or 71 C5 
bulk carriers, and about one million tons of ship's fuel. The same overseas tonnage would 
require over 23 million tons of aircraft fuel and more than 1600 very large aircraft if it were 
to be carried by air. The magnitude of the savings from continued use of water transportation 
is changed little by such refinements as allowing for longer sea routes to avoid Suez. 


VIII. DISCUSSION 


This study has examined the feasibility of a very large scale global airlift in terms of 
fuel and aircraft requirements. The results show that a complete airlift is not practical with 
near-future aircraft and the present demand for strategic imports, even if no attempts are 
made to move the rest of our foreign trade. An airlift supported by surface tankers would be 
feasible but expensive and offers little by comparison with a modern surface fleet. 

Many of the problems associated with protection of sea transport are removed or 
lessened when modern, fast ships are used. The 10-knot "notional" ship, and the fleet of 11- 
knot Liberty ships which the notional ship largely represents, are difficult to maneuver and 
cannot outrun a submarine which can surface and circle a convoy to get in position for an 
attack. The Mariner design, however, is a modern 20-knot cargo ship which is harder for a 
submarine to catch and which also spends less time at sea per trip, an important factor in both 
total ship requirements and needed protective effort. The present Mariner fleet numbers only 
28 (including one AKA and two recently converted to passenger service at the expense of most 
of their cargo capacity) but the design has been tested thoroughly, unlike the hypothetical cargo 
aircraft considered. 
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TABLE Ix 
Number of Dry Cargo Ships to Lift One Million 
* 
Long Tons of Payload Per Year a Lift Distance of p(t!) 





Number of Ships 

Lift Distance 
(miles) 10-knot 7900-ton 

"notional" 





Mariner | C5 bulk carrier 








14.98 11.25 5.18 
16.42 11.82 5.51 
17.86 12.39 5.83 
19.31 12.96 6.16 
20.75 13.54 6.48 
22.20 14.11 6.81 
23.64 14.68 7.13 
25.09 15.25 7.45 
26.54 15.82 7.78 
27.98 16.39 8.10 
29.42 16.96 8.43 
30.87 17.53 8.75 
32.32 18.10 9.08 
33.76 18.67 9.40 
35.20 19.25 9.72 
36.65 19.82 10.05 
38.10 20.39 10.37 
39.54 20.96 10.70 
40.98 21.53 11.02 
42.43 22.10 11.35 














10 20 16 





T 7900 10,000 








*Table IX is computed from the formula 


1,000,000 
T [365 + (2D/24S + 39)] 





where T = cargo deadweight tonnage, D = lift distance, and 
S = sustained speed in knots. 


(1) Althoughthis assumes allcargomovinginthe same direction, 
the requirements also apply for one million tons of cargo 
each way over lift distance D. 





FEASIBILITY OF A GLOBAL AIRLIFT 103 


The full-scale sealift would require 846 Mariners, or 30 times as many as now exist; 
tanker support for the full-scale airlift would need only 261 to 288 T-2 tankers, which could be 
met by the U. S. tanker fleet if mobilized for this task. Construction of 820 Mariners, how- 
ever, would have much more peace-time value than construction of the 5000 very large trans- 
port aircraft. Protection of the Mariner fleet would be much easier than the threefold dif- 
ference in ship requirements (between tankers and Mariners) indicates; much of this additional 
requirement comes from the greater loading and unloading time a cargo ship requires. If 
Mariners and equivalent bulk cargo ships carried all our commerce, the number of ships at 
sea on an average day would be less than twice as many as the number of tanker ships needed 
to support an airlift of the same cargo. This difference is the net saving which creation and: 
maintenance of a dry-cargo airlift would yield in antisubmarine effort required. 

To assist in extending these conclusions, the remainder of this section briefly presents 
some other factors which will influence the utility of an extensive airlift. 


COSTS 

The cost of the airlift analyzed in this study was determined only in terms of numbers 
of aircraft and quantity of aircraft fuel required; the financial costs were not considered. 
Although the monetary costs could be ascertained, they would be largely meaningless, since 
resources rather than money are the limiting factors in a controlled, wartime economy. But 
even measurement of the cost of such an operation in resources is limited in accuracy by 
several elements. 

The two principal limitations result from the effects of rationing in wartime, and the 


use Of low grade domestic resources. At any time a sea lift proved impossible, it might be 
far easier to exploit low grade domestic resources than to import higher grade foreign 


resources by air. Both rationing and extended use of low grade domestic resources appreciably 
reduce the need for imports. But the size of this reduction is impossible to determine, since 
the economic factors governing the use of low grade resources; and domestic rationing cannot 
be estimated. 

As a result, it is impossible to determine meaningful relative costs for large air and 
sea lifts, either in war or peace. For the same reasons, the requirements for numbers of 
aircraft and amounts of aircraft fuel are not entirely accurate, but they do provide a clear 
basis for showing the huge advantages of maintaining and protecting an adequate sea lift. 


IMPORT TRENDS 

The information on imports used in this study applies to the early 1950's, and may, of 
course, be considerably changed within a decade or two. Table XI presents some indication of 
the trends in imports. Although the materials covered in table XI do not correspond exactly 
with those in table III, it covers most of the total tonnage in table III, and includes some other 
vital items such as petroleum. Table XI clearly shows that our requirements for imports, 
particularly for iron ore and petroleum, are growing rapidly. The very large increases in 
imports of heavy materials shown in table XI strongly reinforces the arguments against a 
global airlift, and shows that even with expanding use of airlift, the demand for sea lift will 
continue to increase in the next few decades. Manganese imports have increased even more 
rapidly than the estimate in table X. 
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TABLE X 
Summary of Aircraft and Ships Needed and Relative Requirements 





Cargo Ton-Miles Full Scale Air or Sea Lift | Strategic Materials Lift 





Requirements: 
For complete airlift: aircraft 10,950 to 16,010 1630 to 2290 
For sea-fueled airlift: aircraft 4380 to 5000 736 to 833 

tanker ships 261 to 288 40 to 44 
For complete surface lift: 
10-knot "notional" ships 1390 238 
or C4-S-1a Mariners 846 145 
or C5 bulk carriers 71 


Equivalents: 
Complete vs. sea-fueled airlift: 1 ship = [ 1 ship = 
tanker aircraft vs. ships 25.2 to 38.2 acft. | 22.4 to 33,1 acft. 
Complete airlift vs. sealift: 
10-knot "notional" ships if essed ‘ap- 
| 7.9 to 11.5 acft. 6.8 to 9.6 acft. 
1 Mariner = 
12.9 to 18.9 acft. 


1 Mariner = 
11.2 to 15.8 acft. 


1 C5 ship = 
23.0 to 32.3 acft. 


C4-S-la Mariners 


C5 bulk carriers 
Sea-fueled airlift vs. sealift: 


Sealift by 10-knot ships 32 to 36 acft. + 
1.9 to 2.1 tankers 


31 to 35 acft. + 
1.7 to 1.8 tankers 
10 Mariners = 
Sealift by Mariners 52 to 59 acft. + 
3.1 to 3.4 tankers 


10 Mariners = 
51 to 57 acft. + 
2.8 to 3.0 tankers 
10 C-5 Ships = 
Sealift by C-5 bulk cargo ships 104 to 117 acft. + 
5.6 to 6.2 tankers 


10 ten-knot ships = 10 ten-knot ships = 

















It should be noted however, that the import-export flow and the tonnage-distance "mix" 
used in these computations reflect the 1952 and 1953 points of origin and destination for our 
imports and exports. In the event of war these points of origin and destination and the cor- 


responding distance "mix" should change considerably. No attempt has been made to project 
these shifts. 
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TABLE XI 
Trends in Strategic Imports* 





Percent Increase Quantity Increase 
Expected: 1950-1960 (long tons) 





Copper (crude and semi- 
manufactured) 31.7 175,000 


Lead (crude and semi- 
manufactured) 20.0 103,000 
Zinc (crude and semi- 
manufactured) 30.8 114,000 


Tin (metal and metal 
content of ore) 5.0 5000 


Aluminum 77.6 173,000 


Bauxite 97.5 2,440,000 
Nickel 30.9 27,000 
Manganese ore 28.1 225,000 
Iron ore 300.0 24,000,000 
Industrial Diamonds 50.0 


Other non-ferrous metals 
and ores 19.7 


Crude petroleum 103.0 24,600,000 
Rubber, crude 10.9 87,000 
Residual fuel oil 56.6 9,550,000 
Wood pulp 20.0 434,000 


Newsprint 28.0 1,240,000 
Sawmill products 20.1 

















*Source: America's Needs and Resources, A New Survey, J. Fredrick 
Dewhurst and Associates, Copywrited by The Twentieth Century 
Fund. 1955. In some cases quantity figures are for metallic 
content ofore althoughthis is not specifically stated (manganese 
is an example). These estimates have already proved con- 
servative in some cases (manganese is again an example, with 
an 18 percent increase from 1952 to 1953 alone) but should be 
adequate to show the upward trend. 
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NUCLEAR POWER 

Although the demand for imports is shown to be increasing there is still some question 
as to whether or not innovations in aircraft design and propulsion, such as nuclear power, may 
make aircraft more favorable competitors. The development of nuclear powered aircraft can 
greatly reduce the numbers of aircraft needed for a global airlift, since there will be no need 
for tankers. However, nuclear aircraft will be at least as expensive to build and operate as 
conventional aircraft, and probably a great deal more expensive. 

For the next decade or two, therefore, the best that could be expected from the success- 
ful development of nuclear aircraft would be the situation presented in section VI, which cor- 
responds to an airlift which requires no aircraft tankers, since ships lift the fuel. In this case, 
more than 4000 very large nuclear aircraft would still be needed for a full global airlift. 
Although nuclear power can provide unique advantages to a logistic aircraft, perhaps exceeding 
the advantages to combat aircraft, it seems unlikely to effectively influence the advantages of a 
sea lift in the next two decades. Even in the more distant future, the influence of nuclear air- 
craft power will depend on cost and resource factors which cannot now be estimated. 

In summary, the increasing dependence of this country on the import of heavy materials 
will continue to demand strong dependence on the sea as a commerce carrier, despite tech- 
nological progress in aircraft design and propulsion. In fact, our dependence on the sea is 
certain to increase in the foreseeable future. While airlifts of selected items on a moderate 
scale can handle certain critical materials, meet emergency military needs, reduce require- 
ments for varied inventories at overseas bases and otherwise provide a useful supplement to 
sea transport, this supplemental role is the only one which can be justified for airlift even with 
large cargo aircraft of the types currently proposed. 
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APPENDIX A 
THE DERIVATIONS OF AIRCRAFT FUEL AND SORTIE COST RELATIONSHIPS 


THE AIRCRAFT FUEL COST RELATIONSHIP 

The aircraft and fuel requirements shown in tables I and II in the text are based on the 
performance characteristics of a hypothetical 1963-1968 turboprop super cargo/tanker aircraft 
of somewhat higher performance than any yet designed. The assumed structural character- 
istics of the hypothetical aircraft are as follows: 


Dimensions 
Wingspan: Weight Empty: 176,000 lbs. 
Length: Basic: 180,000 lbs. 


Height: Power Plant: 4-Turbo- Prop 
Improved T-57 
Tread: Type 


Average 
Flight Speed: 400 knots 


For the purposes of this analysis it will be assumed that on a typical mission this air- 
craft will take off weighing 450,000 lbs. and is capable of carrying 120 long tons of cargo 
and/or fuel, or any combination thereof. The average fuel consumption during a typical 
mission is assumed to be at a rate of 0.015 long ton per mile. This is admittedly a simplifica- 
tion which removes many factors of importance, but accuracy to within 10 percent or 20 per- 
cent is adequate to present a quick estimate of the demands a global airlift would make on our 
resources and production facilities. Many refinements would be necessary for computing such 
things as logistic requirements for specific airfields along the various airlift routes. 

The optimum sortie length is affected by several factors which are not considered in 
the 0.015 long ton per mile fuel use rate, such as fuel needed to climb to altitude and additional 
route mileage if airfields are not in a straight line. Under the 0.015 long ton per mile fuel 
consumption rate, short steps are more efficient as they avoid unnecessary transportation of 
fuel. (This shows up in two ways. If fuel must be transported by the airlift, shorter legs mean 
that less fuel is carried one way by tanker aircraft and carried several hundred miles back in 
the aircraft which uses the fuel. If fuel is available without airlift, more frequent fuel stops 
mean more cargo capacity as cargo aircraft fuel loads decrease. The first of these effects 
compounds significantly because intermediate stages of the airlift are transporting less fuel 
and moving it more efficiently.) Table A-I shows fuel consumption per ton of cargo moved 
under three fuel-transport schemes and three step-lengths: 2000, 1000, and 500 miles between 
airports. 

The maximum work (ton-miles of payload) accomplished per aircraft sortie is obtained 
from 2000-mile sorties. For shorter sorties, the work per sortie decreases, but the cost per 
sortie (in fuel consumption and in aircraft hours) also decreases. For longer sorties, fuel 
needs are so great that payload ton-miles per sortie become less even as the cost of a sortie 
becomes greater. A 2000-mile sortie is, therefore, the maximum worth considering, unless 
support airfields cannot be set up at 2000-mile or shorter intervals. The maximum unavoidable 
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TABLE A-I 
Fuel Consumption Per Ton of Cargo Airlifted 


Tons of Fuel Required to Transport One Ton of Cargo 





Full cargo load, Fuel andcargo moved Fuel moved by 
fuel airlifted in opposite directions surface lift 





2000- 1000- 500- 2000- | 1000- | 500- 2000- | 1000- | 500- 
mile mile mile mile mile | mile mile mile | mile 
steps steps steps steps | steps | steps steps | steps | steps 





500 .143 143 143 133 133 133 133 .133 .133 
1000 333 .333 306 286 -286 -267 .286 286 267 
1500 -600 -324 -493 -462 -419 -400 -462 -419 -400 
2000 1.000 778 -706 -667 571 -533 -667 571 -533 


2500 1.286 | 1.032 -950 800 -704 -667 800 -704 667 
3000 1.667 | 1.370} 1.228 953 857 800 953 1857 800 
3500 2.200 | 1.709 | 1.547 1.129 -990 -933 1.129 -990 933 
4000 3.000 | 2.160 | 1.910 1.333 | 1.143 | 1.067 1.333 | 1.143 | 1.067 


4500 3.572 | 2.612 | 2.326 1.476 | 1.286 | 1.210 1.467 | 1.276 | 1.200 
5000 4.333 | 3.214] 2.801 1.667 | 1.476 | 1.373 1.619 | 1.429 | 1.333 
5500 5.400 | 3.817 | 3.344 1.933 | 1.667 | 1.560 1.795 | 1.561 | 1.467 
6000 7.000 | 4.619 | 3.965 2.333 | 1.921 | 1.773 2.000 | 1.714 | 1.600 


6500 8.143 | 5.423 | 4.674 2.619 | 2.175 | 2.017 2.133 | 1.847 | 1.733 
7000 9.667 | 6.492) 5.485 3.000 | 2.513 | 2.295 2.286 | 2.000 | 1.867 
7500 11.800 | 7.563 | 6.411 3.533 | 2.852 | 2.614 2.462 | 2.133 | 2.000 
8000 15.000 | 8.989 | 17.470 4.333 | 3.303 | 2.977 2.667 | 2.286 | 2.133 


8500 17.286 | 10.417 | 8.679 4.905 | 3.612 | 3.393 2.800 | 2.419 | 2.267 
9000 20.333 | 12.319 | 10.062 5.667 | 4.357 | 3.868 2.953 | 2.571 | 2.400 
9500 24.600 | 14.224 | 11.642 6.733 | 4.817 | 4.411 3.129 | 2.704 | 2.533 
10,000 31.000 | 16.758} 13.448 8.333 | 5.762 | 5.032 3.333 | 2.857 | 2.667 






































flight, fortunately, is only slightly over 2000 miles, from California to Hawaii (San Francisco- 
Honolulu equals 2080 n.m.). Costs based on 2000-mile legs are therefore maxima. Full 
savings from shortening individual legs to 500 or even to 1000 miles cannot usually be 
achieved; many inter-island hops are in the 1200-1800 mile range and further reduction 
involves additional mileage to reach islands not on the straight-line route. The only routes 
involving no over-water flights of at least 1200 miles are the South American routes of the 
global airlift. Fuel consumption rates in the body of this study are therefore computed on both 
1000 and 2000 mile legs and should actually lie somewhere in-between. 

Three airlift schemes have been considered. The most costly (in fuel consumption) 
involves airlifting fuel loads in the same direction as full cargo loads. This will apply whether 
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aircraft return lightly loaded or carry capacity loads in both directions. An airlift between 
North America and Europe would be in this category. At the other extreme is an airlift in 
which fuel is available at all airports without airlift. This applies if tanker ships are used, 
and also for limited portions of some likely airlift routes which pass through oil-producing 
areas, or airports which can be served by pipeline or surface transport from an oil-producing 
area. The South American airlifts can be in this category, particularly on the Brazil run via 
Venezuela. An intermediate situation is that in which cargo is moving only one way, in the 
direction opposite to fuel, so that cargo aircraft can be fully loaded with fuel and serve as their 
own tankers. A strategic mineral imports airlift, as discussed in this study, would fall into 
this general category except that the raw material producers would probably demand airlifted 
goods in payment. 

Table A-I is based on the following formulas for these three airlift schemes: 

a. Maximum capacity utilized for cargo; fuel transported by tanker aircraft. Airlift 
involves N steps of equal length K, requiring F tons of fuel for a one-way flight 
over one step. If necessary a shorter step is added of length L. Total aircraft fuel 
required per long ton of payload airlifted in the direction of heavier traffic is 





F (60)272 
=, (60 - F)* 


H + (1 + H) long tons 


where H = the fuel tonnage required to transport one ton of payload over the shorter 
step of length L. 


2 L (0.015) L 
~ 120 - 2L (0.015) 4000 - L 





and is tabulated as follows: 


L 500 1000 1500 
H .143 .333 = .600 


- Maximum capacity utilized for cargo; fuel available at all airfields without airlift. 
Airlift again involves N steps of length K, requiring F tons of fuel for each one- 
way step of K nautical miles. If necessary a shorter step of length L is added. 
Total aircraft fuel required per long ton of payload airlifted in the direction of 
heavier traffic is 





2 NF 2 FL 
120- F * 120K - FL 





. Maximum capacity utilized in one direction and all fuel airlifted in the other 
direction. Airlift again is in N steps of length K (each step requiring F tons of 
fuel each way) plus a fractional step, if necessary, of length L. Total aircraft fuel 
required per long ton of payload airlifted is 
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a=1 (60 - F)* 


with Q= + and an integer less than N, representing the number of steps over 
which no tanker airlift is necessary because cargo lift needs in one direction are at 
L 

4000 - L 
If Q is not an integer, the fuel requirement may be approximated by applying the 
formula twice, using the integers next above and below Q, and interpolating between 
the two resulting fuel use values. If N <Q, then the airlift requires no tanker sup- 
port and the formula of case 2 applies; if N = Q the fuel use value is simply 


least as great as fuel airlift needs inthe other direction. As in case 1, H = 


ae 
120 - F 


AIRCRAFT REQUIRED 

Total aircraft mileage for lifting a given tonnage for a given distance can be obtained 
easily from table A-I, since the table is based on a uniform rate of fuel consumption of .015 
long tons per aircraft mile. Each ton of fuel consumed represents 66.7 aircraft miles. 

If it is possible to maintain an average of 200 flying hours per month and an average 
forward speed of 400 knots, so that the average aircraft produces 80,000 miles per month 
(960,000 miles per year) of airlift mileage, then each ton of fuel consumed per year represents 
1/14,400 of one aircraft-year. Aircraft requirements for any level of airlift can be obtained 
once annual fuel requirements are obtained from table A-I, by dividing annual fuel required, in 
tons, by 14,400. 
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1. INTRODUCTION 

The games examined in this paper were originally formulated by one of the authors in 
order to investigate some problems in economic theory pertaining to the theory of the firm [4]. 
A corporation has to give consideration both to its prospects for survival and to its ability to 
pay dividends. Survival per se could be a goal, as could the maximization of the discounted 
value of the dividend payments. The latter might even involve the eventual ruin of the firm as 
part of the optimal policy. 

In order to portray the relationship between the market and the financial aspects of the 
firm, we can construct a type of game whose form is closely related to the Gambler's Ruin 
problems. The fortunes of the firm are divided into a corporate and a withdrawal account. 
Ruin conditions apply to the former, and dividend payments are made in the latter. 

In general, complex dynamic problems are soon encountered. This paper deals only 
with simple examples involving one and two person dynamic nonzero sum games. 


2. ONE-PERSON GAMES OF ECONOMIC SURVIVAL 
2.1. General Description of One-Person Games of Economic Survival 

A one-person game of economic survival is characterized by the following quantities: 
a discrete set of times t = 0, 1, 2, ...; a set of integers aj, i=1,...,n that represent the 
amounts that a player obtains each time he plays; a set of positive numbers Pi» cs oreo 
with 2 P; = 1, where p; represents the player's probability of obtaining a; when he plays and 
where successive plays of the game are assumed to be independent; a corporate account whose 
size at time t is denoted by C(t); a withdrawal account whose size at time t is denoted by 
W(t); a discount rate p (<1) which is effective on the withdrawal account only; a ruin payment 
F that is paid to the player upon ruin or on leaving the game by withdrawing all assets; a ruin 
level B and a ruin condition that says, if C(t - 1) > B and C(t) < B the player received the 
ruin payment F and the game is terminated at time t. 

At the beginning of any time period he may transfer any integral part of his corporate 
account into his withdrawal account; let us indicate such a transfer by w. Such a payment may 
be interpreted as a dividend payment to stockholders and is nonrecoverable, hence payments 





*Manuscript received September 12, 1958. 

+The authors wish to express their gratitude to Drs. Lloyd Shapley and Alan Hoffman for 
valuable criticism. Additionally, the authors are indebted to Mr. R. Singleton for several 
interesting comments and observations, including an alternative proof for the solution to the 
one-person game. 
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can never be made from the withdrawal to the corporate account. Let the initial amount in the 
corporate account be C(0) = x (an integer). The goal of the player is to maximize the dis- 
counted expected value of his withdrawal account. This completes the description of the game. 

Suppose for a moment there were no transfer privileges. Then if we watch the corpo- 
rate account it varies as if it were a random walk on the line. Thus it moves from a value 
C(t) at time t to the value C(t + 1) = C(t) + a; at time t + 1 with probability p;. 

If there are transfer privileges, then, since we have assumed that each time the game 
is played the probabilities of obtaining payments a; are independent of previous plays of the 
game, we see that the transfer payment w should be a function only of the value of the corpo- 
rate account. Therefore, if transfer payments w(C(t)) are included, the corporate account 
moves from C(t) to C(t + 1) = C(t) + a; - w(C(t)) with probability p; for i=1,°**,n. 

A financial strategy (or simply, strategy) is a transfer payment function w(x). (The 
function w is from positive to non-negative integers.) Let W be the set of all financial 
strategies. Given a strategy w, one can compute the value function V(x, w) that measures 
the worth to the player of playing the game with initial assets x and strategy w. The value 
function gives the expected discounted value of the player's withdrawal account. 

In order to compute the expected discounted value of the withdrawal account, given 
initial assets x and withdrawal strategy w, we consider truncated games of length n, that is, 
games which are played at times t=0,1,°°*,n. Let Sh be the outcome of a series of plays 
of the subgames played each period. This indicates whether the player has obtained the maxi- 
mum or minimum revenue from the first play of the game, the second, etc., up to the nth play. 
From this information we can compute the amount K(s,, w,t) that is added to the player's 
withdrawal account for each of the times t=0,1,°**,n. Let q(s,) be the probability that S, 
occurs. Then, taking sums and passing to the limit, we have 


n> Ss 
as n 


n 
V(x,w) = lim )) pa K(s,,w,t) p* q(s,) 


\t=0 


as the expected discounted value of the withdrawal account. In order to show that this limit 
exists, we consider the most favorable game—the one in which the player wins every time he 
plays. For this game, the most favorable strategy is w, for which the player withdraws all 
his assets above the ruin level as soon as he can. Then 


[°.¢) 
V(x,w') = (x - B) ‘i pf = (x- B) +P 


for this game. And, obviously, for any other course of the game and any other strategy w, we 
have 0 < V(x, w) < V(x,w ). Since the expression for V(x,w) is a series of positive terms 
and is bounded, it converges. For an explicit calculation of V(x, w) in the case that w is an 
optimal strategy, we refer the reader to Section 2.3. 

A strategy w is said to dominate strategy w if 


(a) V(x,w)>V(x,w) forall x>1; 


(b) V(x,w) > V(x,w') for some x>1. 
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That undominated strategies exist can be seen from the following argument. A with- 
drawal strategy w is completely characterized by the smallest integer N such that w(N) > 0 
since the corporate account increases only by integral values and N acts as a "reflecting 
barrier." We shall call this N the N-value of w. Now we have K(s,,, W, t) = 0 for 
n=1, 2,...,N-x, so that, letting ¥,. (x, w) be the discounted amount in the withdrawal 
account out to time n, we have 


V, (% W) = 2 i K (s,,, w, t) p' a(s,) 
=0 


n 


- {= K(s,, w, t) q(s oe 
=0 


n 
en 7 z Kley 49 a6) pi Hn 


t=n-x n 


Now the term in braces on the right is bounded by the same kind of argument that we used to 
show that V(x,w) was bounded. Since p < 1, the factor pN-x becomes small if x is fixed and 
N becomes large. Hence all those strategies w whose N-values are sufficiently large are 
dominated by other strategies with smaller N-values. Since there are only a finite number of 
withdrawal strategies with bounded N-values, there are only a finite number of values V(x, w) 
to consider, and hence an undominated strategy w exists. 

Let W° be the set of undominated strategies. In Section 2.3. we shall characterize w 
for a certain kind of one-person game of economic survival. 


2.2. A Simple One-Person Survival Game 

Consider the game G for which a= i, ao = - 1, Py? 0, Po > 0, Py +Po= 1, B=0, 
and F = 0. We shall characterize the undominated strategies for this game and thus obtain 
information about the result for the general one-person survival game. 


LEMMA 1: An undominated w® in G is a solution of the functional equation* 


(1) V (x, w°) = max [ Py V(x - 1, w°) + PP» V(x +1, wy), max [a +V(s-a, wy] . 
a>l 


*The optimal financial strategies are not the only solutions of (1). For instance, the function 


x 
1+//1 - 492 
vee>| ae 


2Ppp) 





satisfies the equation 


Vid = ons Pp, Ve- s+ 0+ Pp, Ve ~ 0+ Ute). 
a>0 


and we can regard this function asthe solution for a zero withdrawal strategy. But this solution 


is not economically interesting, since, although V (x)—>o, the expected discounted value of the 
withdrawal account for any finite length of time is zero. 
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PROOF: . Equation (1) states merely that the player always has the choice of imme- 
diately playing the game or first withdrawing an amount a and then playing. Moreover, 
he must choose the withdrawal payment a so that his value is maximized. It is obvious 
that any strategy, not a solution of (1), is dominated by such a solution. 


COROLLARY: V[x, w’] is a monotone strictly increasing function of x. 
PROOF: From (1) we have V(x, w°) > 1+ V(x - 1, w®), hence V(x, w°) > V(x - 1, w°). 


Lemma 1 shows that the game is a dynamic programming problem in the sense of 
Bellman [1]. It also brings out a difficulty: namely, suppose that for one or more values of a, 
both terms of the right-hand side of (1) are equal. An undominated strategy could then choose 
either term. We determine a unique choice by saying that he will not play the game when the 
terms are equal. More precisely we give the following definition. 


DEFINITION: By a special undominated strategy we shall mean an undominated 
strategy with the following property: wo (x) = 0 if, and only if, 





(2) PP, V(x - 1, w°) + PPy V(x +1, w’) > max [a + V(x - a, w°)] . 
a>l 


Remark: If p Py < 1/2, the solution to the game is trivial, the player pays out all the resources 
to his withdrawal account immediately as the game is "unfair" and not in his favor. We exclude 
this case and henceforth assume Pp Py 2 1/2. 


LEMMA 2: If w® is undominated, then there exists an x such that w? (x) > 0. 
PROOF: If w° (x) = 0 for all x, then, obviously, V(x, w’) = 0. But, since all assets 


can immediately be withdrawn, we see that V(x, w’) >x >0, which is a contradiction 
for positive x. Therefore for some x, w? (x) > 0. 


LEMMA 3: Let w® be a special undominated strategy. If z-is an x such that w? (x)=0 
and w°(x +1) > 0, then w°(z +1) =1. 


PROOF: Suppose, on the contrary, that w? (z+1)=k+1>1. By Lemma 1 we have 
V(z +1, w’) > V(z,w°) +1 >V(z - 1, w°9) +2 >... >V(z-k, w) +k. 


Since the player always withdraws to produce the maximum amount possible, by Lemma 1 we 
have 


V(z +1, w))=V(z+1-n, w) en, 
hence all of the above inequalities are actually equalities. And, in particular, V (z, w°) 


= V(z - 1, w°) +1. The definition of a special undominated strategy now shows that w°®(z) > 0, 
contradicting the fact that w°(z) = 0. 
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LEMMA 4: Let w® be a special undominated strategy and let z be the smallest x 
such that w® (x + 1) > 0, then for all x >z we have w? (x) =x- Zz. 


PROOF: No payment can be greater than x - z, for if w? (x)=a>x-zthenz>x-a 
and 


V(x -a, w’) +a >V(z, w°) +(x - z), 
or 
V(x - a, w°) + (z - (x - a)) = V(z - (z - (x - a)), w°) + (z - (x - a)) > V(z, w®), 


contradicting the fact that w°? isa special undominated strategy and z is the smallest 
x such that w(x + 1) > 0. (See (1).) 


If there is at least one x >z such that w® (x) <x - z, then there must be an x > z such that 
w° (x') = 0, for otherwise (2) would be violated. (For instance, x' = x - w? (x) has this 
property.) Let z be the least x > z such that w(x) = 0. By an argument similar to that of 
Lemma 2, there is an x > Z such that w(x) = 0 and w(x +1) > 0. Call the smallest such Z4- 


By an argument analogous to Lemma 3, w(z, +1)=1. Then we have the following information 


about w’: 


0 for x <z 


x-Z forz<x< Z 
Oo 
w (x) = o 

0 forz <x < Zy 


1 forx=2,+1. 


We shall show that w°® is actually dominated by another strategy w' defined as follows: 


(f QO forx<zor ifx= 2, for the first time, 
x-z for (i)z<x<z-1 if x has always been < z, , 
(ii) z<x<Z if x has ever been =Z,; 
for (i)Z-1<x<z, if x has always been < z, , 


(ii) Z <x £24 ifxhasever been =2Z,. 





L wo (x) forx>Z,+1. 


We shall show that strategy w° does not do as well as w’ by comparing two games, 
one starting with C(o) = a and using strategy w’, and the other starting with C(o) = a +1 and 
using strategy w’, where Z <a<a+l1< Z1: Let C’ (t) be the player's corporate account in 
the first game and Cc (t) be his corporate account in the second game. We compare these two 
strategies for the same course of outcomes of individual plays of the game. 








M. SHUBIK AND G. L. THOMPSON 


For some t >1, C' reaches Zy before C° reaches Z - 1. Then 1 is with- 
drawn from the C° account, but the C' remains unchanged “ that at time t 
the two ee accounts will be equal to Zy- And since w' becomes the 
same as w° after reaching Z1> the two corporate accounts will be tag for 
the rest of the game. Hence we see that V(a +1, w °) - V(a, w ‘y=p' <1. 


Case 2: For some t > 1, C° reaches Z - 1 before C' reaches zy. Then Z-1-z is 
withdrawn from C° and Z - 2 - z is withdrawn from C' so that at time t the 
two corporate accounts will be equal (to z); and since w' is the same as w® 
when x < Zz, the two corporate accounts will remain equal for the rest of the 
game. Hence we see that V(a + 1, w’) - V(a, w) = pt <1, 


In both of the above cases it is false that V(a + 1, w°) > V(a, w’) + 1, which must 
necessarily be the case for an undominated strategy. (One could do better than w? at a+1 by 
first paying out 1 and then using w'.) From this contradiction we conclude that Z does not 
exist and w°(x) = x - z for x > z, concluding the proof. 


THEOREM 1: There is a unique special undominated strategy w° defined as follows: 
there exists a z >0 such that w° (x) = 0 if x <z and w° (x) =x-zif x>z. 


PROOF: Observe from the definition of special undominated strategies that z (the 
largest x for which w? (x) = 0) is the same for all special undominated strategies. Lemmas 
3 and 4 and mathematical induction now show that any two special undominated strategies are 
equal. 


Remarks: 

(A) If p= 1, then this one-person game becomes the classical gambler's ruin game. 
The unique special undominated strategy for this game would be to initially withdraw all funds 
and not play the game at all. 

(B) The calculation of the value of an example of a game of the type here discussed is 
carried out in Section 2.3. 

(C) If the changes of state are rational numbers ay and ao, where a; > 0 and ao < 0, 
then the changes in the possible fortunes of the player will be multiples cf d, where d is the 
greatest common divisor of ay and ao. We regard the player as ruined whenever his corporate 
account is less than |@g|. The optimal strategy for the player is as follows: (a) initially pay 
out an amount Y = x - m |ao| where m is the shige integer such that y is mae 
(b) there is a number z such that, if x < z, then w °(x) = 0, and if x > z, then w°(x) = x - z. 
(It can be shown in this case that x - z is a multiple of d.) 


2.3 The Value of a Simple One-Person Game of Economic Survival 

In Section 2.2 we have established the form of the special optimal strategy w for a 
game with ay=+ 1 and ag = - 1, and F=0. The value of a game in which the special optimal 
strategy is used can be calculated directly as follows: Let z +1 be the first value of the cor- 
porate account at which a payment is made to the withdrawal account under the optimal 


strategy w’, By setting Pi =p and Po = q=1- pp, the value of the game is the solution to: 













ca 
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(3) V (x, w°) =ppV(x +1, w’) +pqV(x- 1, w’) for x<z 
subject to: 

(4) V(0, w°) = 0, 

(5) V(z +1, w’) = V(z, wy) +1. 


x 


We set V(x, w’) = a” in (1), which yields: 


i: 1+V/1- 4p?pq_ 


7 2pp 





Denoting the roots by ay and Qo, 
V (x, w’) = Ay at + Ay as " 
From (4) and (5) we obtain: 


Ay 


- Ao 
and 
1 


z+1 Z+1 z z\ 
(a3 “oS ) “ (a3 - af | 


The solution may be expressed as: 





A, = 


V(x, w°) = A; (z) for - a3} ‘ 


The solutions for games with pp > 1/2 for different values of p and p have been 
calculated and examples are provided in the tables below. 
































p=.7T, p= .99 p= .995, p= .99 
P 
Zz Zz 
x x 
1 2 3 2 3 4 

1 1.59 1.10 1 191.34 195.08 195.08 
2 2.29 1.60 2 194.21 198.02 198.02 
3 2.29 3 199.04 | 199.04 
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p = .99, p= .55 p = .995, p= .55 














10 








3.06 
5.99 
7.71 
9.52 
11.09 
12.47 
13.72 
14.87 
15.95 


The numbers in the boxes represent the values of the largest investment that will be 


made under an optimal policy. The underlined z's are the amounts associated with the optimal 
policies. 
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3.1. GENERAL DESCRIPTION OF TWO-PERSON GAMES OF ECONOMIC SURVIVAL 

A two-person game of economic survival is characterized by the following quantities: 
an mxn matrix A= lla. ||, a discount rate p (<1), ruin payments Fi and Fo, ruin conditions 
B, and Bo, corporate accounts C = (C, (t), Co (t)), and withdrawal accounts W = (W, (t), Wo (t)). 
The role that each of these quantities plays is analogous to those played by the corresponding 
quantities in the one-person case. The subscript 1 refers to the first player and 2 refers to 
the second player. Let the initial amounts in the players’ corporate accounts be C, (0) =x 
and Cy (0) = y. We shall seek equilibrium point solutions for a simple class of these games 
for various values of x and y. 


3.2. A Simple Two-Person Survival Game 
Here we consider the case in which B, = Bo = 0, Fi = Fo = F, and A is the matrix 


We shall characterize the equilibrium point solutions of this game for various initial fortunes 
xX and y. 

In general it is false to assume that the players will always use optimal strategies for 
every matrix game A, as will be shown in Section 3.4. In the simple example selected here, 
however, the unique optimal strategies for each player are to choose each of his alternatives 
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with probability 1/2, and by known martingale theor 2ms [2] there is no optimal way to deviate 
from these. Hence each player can expect +1 and -1 with equal probability each time the game 
is played. 

Let us compute the values V(x, y) = (V, (x, y); Vo (x, y)) of withdrawal accounts of each 
player, assuming that they begin with fortunes x and y and that they both always keep their 
fortunes in their corporate accounts. 

We must distinguish between two cases: the first, where at the termination of play the 
winner obtains only the "prize'' F, and the second, where he obtains the prize F and adds the 
assets in his corporate account to his winnings to bring his total to x + y + F. Initially, we 
examine the first case. However, it is easier to give an economic interpretation to the second. 
This is done later in Section 3.4. 

For the first player the quantity Vi (x, y) satisfies 


(6) V(x, y) =2Vy (x +1, ¥- 1) +E Vy (x- 1 y+). 


By use of the trial solution V, (x, y) = c*dY where c and d are constants, and the boundary 


conditions V, (x + y, 0) = F and Vo (0, x + y) = 0, it can be shown that the solution to these 
difference equations is 


x x 
7 V F ‘1 om 
(7) 1% y) = X+y x+y 5 
~ i 


where 


1-1 - p* 
(8) me awn aaibbintiinicde , = ———_= 


p 


i 
+ 


If 0<p <1, then c, > 1 and 0<cy <1. (This is certainly the range of economic interest.) 
If x +y is relatively large, then the second factor in the denominator of (7) is negligible and 
we may write 


F "9 
(9) Vv; (x, y) = —(1-k*) where k =— = ce (k <1). 
c} "5 


We use this approximate relation to investigate the solutions. Note that for fixed y the 
quantity Vv, (x, y) is a monotone increasing, but bounded, function. Thus, given a fixed fortune 
y for player 2, player 1 will not want his capital account to be arbitrarily large, since it would 
be more profitable for him to withdraw some of this money and invest it elsewhere. 

Now suppose the players can withdraw any part of their corporate account, so long as 
neither player is ruined. But if either player is ruined, the survivor gets F and forfeits the 
money in his corporate account. 

Player 1 will withdraw a dollar whenever his expected value obtained by leaving it in is 
less than a dollar more than his expected value obtained by withdrawing it—in other words, 
when (for fixed y) x is such that 
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F F 
(10) ——(i -k + 1 > —-. 
e} e} 


After some algebraic manipulations we take logarithms and obtain a relation of the form 


(11) 


where b is a positive constant. 


Since we have made an approximation in the course of deriving the very simple formula 
(11), we must expect our results from now on io be quantitatively inaccurate. However, for 
large discount rates (small discount factors) the error due to the approximations is small, so 
that our conclusions are qualitatively correct. 

Carrying out the analogous procedure for player 2 we find that, given a fortune x of 
player 1, the maximum fortune y that player 2 would keep in the game is given by the equation 


(12) 


where the constant b is the same as in equation (11). 
We are now in a position to analyze the equilibrium points in the game. 


3.3. Equilibrium Points in the Two-Person Game 


t------- 








Figure 2 


In Figure 1 we have sketched 
the graphs of equations (11) and (12) 
and have numbered different areas of 
the plane that lead to essentially dif- 
ferent behavior, described by a pair 
of strategies in equilibrium. 

In area 1 neither player has more money that he would like to have, based on the most 
pessimistic assumption that his opponent will not withdraw money from the game. It is impos- 
sible for them to move to a point at which either player would prefer to withdraw money rather 
than keep it in. Hence, if the initial fortunes of the players specify a point A in the interior of 
1, as the game progresses their fortunes move along a negative 45° line (see Figure 2) until 








elo R—-—--—-4+--—-—--— ------------ 


i 
y = b-ax 


Figure l 
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the game terminates at one of the axes at which point one of the players is ruined. The equi- 
librium points also have maxmin properties in area 1. 

If the initial fortunes specify a point B in the interior of area 2, then again neither 
player wishes to withdraw money. However, with probability 1, the play of the game will 
proceed until one of the players makes a withdrawal. This is shown in Figure 2 by the 
crossing of the negative 45° line with one of the lines of equation (11) or (12). After one of the 
players withdraws part of his fortune, the point representing the joint fortunes in their cor- 
porate accounts moves back into area 2, and the game progresses until it eventually terminates 
at (0, b) or (b, 0). A possible path of play is illustrated in Figure 2. 

Let A be an initial fortune point in area 3 (see Figure 3). Here the second player has 
more money than he needs under a minmax assumption concerning his opponent's behavior. 

He withdraws down to the point A' which is in area 2. The game then proceeds as if it started 
at A' in area 2. 

In Figure 3, point B is in area 4 where both players are overcapitalized. Two possible 
equilibrium paths are as follows: first, player 2 may withdraw capital so that the fortune move 
to point B' which is in area 2 and the game proceeds from there; or player 1 may withdraw to 
point B'', after which player 2 withdraws to B'"' in area 2 and the game proceeds from there. 
Other equilibrium paths between B' and B'"' are also possible. 

Point C in Figure 3 represents a point in which player 2 has enough capital to make it 
rational for player 1 to withdraw from the game on a minmax assumption concerning player 
2's behavior. This is shown at C''. If, however, player 1 retained his capital, player 2 would 
withdraw down to C'. There is mixed strategy equilibrium involving the decisions of each 
player to keep in ell his funds and the decision of player 1 to leave the game and player 2 to 
reduce investment. There are also many other equilibrium points. 

The equilibrium point solutions in areas 3', 4', and 6' are analogous to those in areas 
3, 4, and 6, but with the roles of the two players interchanged. 

Point A in area 5 of Figure 4 represents a point at which both players are overcapi- 
talized. There are two obvious equilibrium points: in one player 2 withdraws to A' and then 
player 1 withdraws to A'''; in the other, player 1 withdraws to A'' and then player 2 withdraws 
to A'''', There are many other equilibrium points. In particular, E is enforceable as follows: 

player 2 announces that he will withdraw to E' 
and, even though overcapitalized, will remain 
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there, regardless of what player 1 does; it is then rational for player 1 to withdraw to E, and 
the game proceeds from there. There is an analogous equilibrium point in which 1 withdraws 
to E'' and then 2 withdraws to E. 

At point B in area 7 of Figure 4, player 2 is so overcapitalized that player 1 should quit 
entirely. If he does, the fortune point moves to B''. On the other hand, if player 2 first with- 
draws his capital, the fortune point moves to B'; then player 1 should withdraw to B''' and the 
game proceeds in area 2; there is a mixed strategy equilibrium involving these choices. The 
equilibrium point E is also enforceable in the same manner as stated in the paragraph above. 
Many other equilibrium points are possible. The behavior in area 7' is analogous to that in 7, 
with the roles of the players interchanged. 


In area 8 both players are heavily overcapitalized and there are many equilibrium 
points, including E. 


The analysis is most satisfactory in areas 1, 2, 3, and 3', since unique equilibrium 
points are obtained and these have maxmin properties. In other areas, the multiplicity of 
equilibrium points means that the outcome of the game will be settled by means of threats and 
counterthreats, in which case the equilibrium point analysis is not particularly conclusive. 
Experimental investigation of the course of play in areas 5 to 8 would be of interest. 


3.4. An Economic Interpretation of the One- and Two-Person Games 

The one-person game serves to demonstrate the "safety value" of liquid assets in a 
fluctuating market. It is analogous to an inventory problem. The penalty of being "out of 
stock" is the ruin of the firm. 

The next step in the investigation of models with greater realism involves having the 
rewards from the random walk depend upon the asset position of the firm. In the model in 
Section 2.2, the only use of the money in the corporate account was as protection against ruin. 

It is possible to give a mathematical formulation for the policy of a firm which wishes 
to achieve and maintain a given level of safety prior to paying dividends. A discussion and 
examination of the effects of different policies is given elsewhere [4]. 

The two-person game has an economic analogue a market in which there is overcapacity 
with two firms present, but in which one firm could make a profit if it were the sole survivor. 
The zero-sum matrix indicates that if one firm makes a profit, the other must lose. For 
example, the only way in which business may be obtained is by obtaining the competitor's 
accounts. The prize, F, is the discounted value of the income obtainable by the surviving firm 
in the market. When one firm is ruined the value of the other is F plus its corporate assets 
at that time. In order to express this we would have had to express: 


V(x +y,0)=x+yH+F. 


If x + y <<< F, then the same qualitative results as above hold, with the linear equations (11) 
and (12) being replaced by nonlinear equations whose graphs are convex upward curves similar 
to the straight lines of Figure 1. 

If cooperation and side payments were permitted, the solution would be for the industry 
to "rationalize" immediately, with one firm exiting but being paid off by the other. 


An examination of the solution shows that the richer both firms are, the worse off they 
may be!!! 
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V(x, x) > V(x + k, x +k). 


This phenomenon is caused by the prolongation of the length of the struggle for survival. If 
two financially weak firms compete, one will soon fail, leaving the other a lucrative market. 
If they are stronger, it will take a much longer time to reach a prosperous state, unless there 
is collusion and rationalization. 


3.4. Comments and Further Problems 
It is evident that the investigation of the game with: 


-1 
holds for A = 


by a simple change of units. 

If the matrix A is 2x 2 and the one-period zero-sum game on A has a unique mixed 
strategy solution, then the players will always use optimal strategies in the matrix game each 
period. For A larger than 2~x 2, this is not necessarily so. A simple example demonstrates 
this: 


-1 
1 
- 10 
10 - € 


The strategy = . 0, 0) dominates (0, 0,5, 2 in the one period game but not neces- 


sarily in the dynamic game in which the duration of play as well as the expected gain per period 
is of importance. 

A further discussion and economic interpretation of this type of game including cases 
where the one-period "subgame" is nonzero sum is given elsewhere [4]. 
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It is shown that there are precisely nM-lmn-! basesin atransportation 
problem having m destinations and n origins. This number gives an 
upper limit to the number of iterations in solving the problem and 
is much smaller than that given by the combinatorial formula 
(mn)! /(m+n-1)! (mn+l-m-n)! 











It is of interest in any linear-programming problem to know an upper bound N to the 
number of iterations which may be required in the application of the simplex method before an 
optimal solution is obtained. For a problem involving s unknowns and t constraints: 


RE* Ze 


the maximum number of bases is: 


s! 


(1) ~ tl(s-t)! : 


if we assume the rank of A is t. The actual number of feasible bases can be much smaller 
than N, because not all bases need exist or be feasible. 

For a problem of the transportation type, it is interesting to note that the exact number 
of bases can be found. This number is much smaller than that computed from (1). It is the 
very special structure of the matrix of coefficients which makes this computation of the exact 
number of bases possible. 

The transportation problem can be formulated as follows: 
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co ere ee 


minimize : : cj; x 


— ie 


The matrix of the coefficients in the constraints can be written: 


Ih 


where identity matrix of order n, and 3. = (1, 1, 1, ..., 1) n component row vector. 


One can immediately verify that the rank of this matrix is less than (m+n), since 
i first m rows = > last n rows. 


In fact, the rank is precisely (m+n-1), for if we form the determinant D omitting row 
(m+1) and take columns 1, (n+1), (2n+1), ... (mn-n+1), 2, 3, 4, ...n, we have: 








ee pe ee 


In order to find the number of bases in A (i.e., the number of different combinations of 
(n+m-1) linearly independent column vectors), we shall first determine the number of 
determinants of order (m+n-1) in A which are different from zero. We shall make use of a 
theorem of matrix algebra [1, p. 85] which says that: Given a pxg matrix R anda gxp matrix 
S, where p < g, the value of the determinant of the matrix product RS is equal to & Ri Si, R, 
being any pxp determinant from R composed of the columns ij, io, wieg ib and 8; being a 
pxp determinant from s composed of rows ij, reer ib ‘ 


If A is the (m+n-1)x mn matrix obtained by deleting from A the jm row, any 
determinant Dp} of order i (1 < i < m+n-1) extracted from A is equal to 0 or +1. Theproof 


is immediate. Any column of D} contains two 1's, one 1, or no 1's. Let A be partitioned into 
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two sets U, composed of the first m rows, and L, composed of the last n rows. If there are 
two 1's in every column of Di, the sum of the rows of D) which were originally in U, minus 
the sum of the rows which comes from L, gives a row of zeros, and therefore D! = 0. When 
there is no 1 in at least one column of DI, this determinant is also null. When there is a 
single 1 in at least one column, we expand D) by this column, and D) =+ Os Ds being a 
determinant of order i-1 from Al, We can now apply the same reasoning all over again to 
Di. Finally, D} = 0,11, and thus D) = 0 or +1 (for all i, j, 1 <i< m+n-1,1<j < m+n). 
Using the above theorem, we immediately see that the value of the determinant 


: ae : 
(A?) (A) , which is the product of A) and its transpose, is the number of nonzero determinants 


of order (m+n-1) in A), This follows because: 


(5) a) = (a}) (a4) a a a 
a 


2 


the sum being taken over all combinations of the mn columns in A) taken (m+n-1) at atime 
(while maintaining the column order). 


Let us now compute the Al's, j=1, ..., (m+n). Those which correspond to 


j = (m+1),... , (m+n) (that is, when A} is obtained by deleting one of the last n rows in A) 


are identical. Let a) be this determinant. It is: 


E — is an hxk matrix, of which all elements are unity. 


By suitable addition or subtraction of rows and columns, a} can be rearranged as 
follows: 


“ni m-1 


eee ee ee — 














alae ae 
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a} a ypm-1 n-1 ‘ 


The m, A J's which correspond to j=1,...,m also are identical, and the common 


determinant A) has the same value as a}, The calculations are similar to those given above 


and will not be repeated here. 


j 
An = 


Inasmuch as there are n ways to obtain a} and m ways to obtain 7 the total 


number of (m+n - 1)st-order determinants in A is 
(9) (m+n) n™M-1 n-1 | 


Finally, we shall show that, if in an (m+n) x (m+n-1) matrix A formed from 
(m+n-1) columns of A there is one nonzero determinant of order m+n-1, all (m +n) 
determinants of order m+n-1 are different from zero. We have seen that for A any row of 
L (or U) is equal to the sum of all the rows of U (or L) minus the sum of the remaining 
rows of L (or U). This property remains valid for any a also. 

Therefore, if D is a determinant of order (m+n-1) different from zero in A, the 
jth row of A which does not appear in D* is a linear combination of all the rows of D’, 
with each coefficient being either +1. Let us assume that a determinant of order (m+n- 1) 
extracted from ry which includes row i vanishes. This means that a linear combination L=0, 
with not all coefficients vanishing, exists between row i and a set S of (m+n- 2) rows in D’. 
If the coefficient of row i in L is zero, there is a linear relation between (m+n - 2) rows of 
D’ and therefore D* = 0. If the coefficient of row i is different from zero, we replace this 
row by its representation \ as a linear combination of the (m+n-1) rows of D . This gives 
a linear depencence between the (m+n-1) columns of D* with at least one coefficient, the 
row of D* not in S, not zero. Thus, again D =0. In both cases our conclusion contradicts 
the hypothesis D* + 0. Therefore, any determinant of order (m+n-1) from rx including 
row i, does not vanish. 

Inasmuch as the above argument holds for any i, we conclude that if one determinant 
of order (m+n-1) in A is from zero, all (m+n) determinants of order (m+n-1) in A are 
nonvanishing. 

To find the number of bases in A, we must then divide (9) by (m+n); i.e., the number 
of bases is: 


(10) 


The above equation gives the exact number of bases. The number of feasible bases 
will be less than or equal to this number. The aj, b; determine whether any given basis yields 
a feasible solution. 
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AN INTEGER LINEAR- PROGRAMMING MODEL FOR MACHINE SCHEDULING* 


Harvey M. Wagner 


Stanford Universityt 


I, INTRODUCTION 

Several authors [1, 2,3, 10, 12,13] have previously considered versions of the following 
machine-scheduling or sequencing problem: Given n items, each to be processed on one or 
more of m machines, the order of processing for an item being partially or entirely specified, 
find the sequencing of items on the machines which minimizes the total elapsed time to com- 
plete the manufacture of all items. It is assumed that the manufacturing time of an item ona 
machine is specified (i.e., nonstochastic) and that in-process inventory is allowable. 

In this paper we offer an integer linear-programming model which is capable of 
characterizing the problem. We feel it is of interest to present the model, because (1) it 
shows that there is a single model which encompasses a wide variety of machine-scheduling 
situations, (2) the linear format establishes the existence of a finite algorithm which mono- 
tonically seeks an optimum solution to such sequencing problems, (3) although the model is 
now of very limited computational interest, future developments in integer linear programming 
and methods for efficiently handling "secondary constraints" may make numerical solutions of 
particular problems possible, and (4) there is an indication of the possibility of constructing 
special algorithms to exploit the structure of certain of the "classical" scheduling problems. 
Needless to say, the major justification for considering such an approach is that Gomory [9] 
and others [5,6, 11] have recently discovered promising methods for solving integer linear- 
programming problems; a related justification is that Dantzig, Fulkerson, and Johnson [8] have 
achieved noteworthy success in using secondary constraint techniques for solving a problem 
which a priori contains a mammoth number of restrictions. As will be evident below, the 
model in its present form is computationally unwieldy except perhaps for situations with a 
very few machines and a limited number of items; in such cases, a frequently recurring 
sequencing problem or one involving a considerable financial sum might profitably be solved 
by the method herein. We are more hopeful, as stated above, that this presentation might serve 
as a link to a more computationally feasible formulation. 

In Section II we discuss a general model which is adaptable to a large number of 
scheduling situations. The familiar problem [3, 10,12] of sequencing n items, each of which 
must be placed on all m machines, the order of manufacturing being machine 1, machine 
2,..., machine m, is considered in Section III. The further restriction, to three machines 
only, is explored in Section IV. 





*Manuscript received November 20, 1958. 
+ Work reported here was sponsored by the Office of Naval Research. 
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Il. GENERAL FORMULATION 

We may picture the restrictions characterizing a specific scheduling sequence problem 
by means of two tabular arrays. In the first matrix (I), a row corresponds to one of the n 
items, and a column corresponds to a process stage in the manufacture of the item. For 
example, if item i must be processed on machine i, 2,..., m in that order, we define for the 
i-th row the entry at process stage 1 to be "machine 1," at process stage 2 to be "machine 
2," ..., at process stage m to be "machine m." If an item is not to be placed on every 
machine, then the number of process stages is less than m. For the sake of expositional 
simplicity, we postulate that the processing specifications for any item are such that any 
machine does not appear more than once (if at all) in the i-th row of matrix I; a minor 
modification of the model may be made to allow for multiple processing. If at some point of 
the manufacturing process of an item, the item is to be placed on q machines out of a subset 
(Q) of machines, the specific order of manufacturing being unimportant, then we define the 
entry at that process stage to be the subset (Q) and the specification q; e.g., item i at the 
fifth process stage may be required to be placed on any of three (= q) machines in the group 
(Q) = (machine 5, machine 8, machine 26, machine 27, machine 35). Note that all items do not 
necessarily have the same number of process stages, and within each stage we permit several 
of the aforesaid restrictions. Thus we shall consider those production situations in which it is 
possible to analyze the manufacturing process for each item as a consecutive sequence of 
stages, each comprising processing on one or more machines, and such that the order of 
manufacturing within a process stage is irrelevant. By a straightforward extension of our 
analysis, we could also treat situations requiring processing "substages"; to keep the present 
exposition within bounds, we omit the details required for this extension. 

In the second matrix (II) the rows correspond to the items and the columns to the m 
machines. The entry Vs at the intersection of the i-th row and k-th column represents the 
manufacturing time (set-up plus production)* for item i on machine k. If matrix I indicates 
that under no circumstances is an item to be placed on a particular machine, then the cor- 
responding element in matrix II will not be defined. We postulate that our time units are 
measured such that every 1) is an integer. 

We make the following definitions of nonnegative integer-valued variables 


K) J 1 if item i is scheduled in order-position j on machine k 


4 


0 if item i is not scheduled in order-position j on machine k 


n{*) = time at which the item scheduled in order-position j begins processing on 
machine k 


aX) = time elapsing on machine k between completion of processing the item scheduled 


in order-position j and initiation of processing the item scheduled in order- 
position j+1 


ol) = time elapsing on machine k between the initiation of production and the start of 
processing the item scheduled in the first position. 





*We might also let the processing time be a function of the number of items previously 
placed on machine k; we forego this generality as it would further complicate the notation, but 
the interested reader may supply the necessary amendments. 
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By "the item scheduled in order-position j ... on machine k" we mean that, prior to 
this item being placed on machine k, (j- 1) items have previously been processed thereon. 

We let n(k) be the maximum number of items that might ever be processed on machine 
k, This number may be found by scanning matrix I and counting the total number of times 
machine k appears in the process-stage listing. It represents an upper bound to the number 
of ordered positions that need be considered for machine k. We also define N(k) to be the set 
of items for which machine k appears at some stage. 

The first collection of constraints ensures that each item i completes the necessary 
operations within every process stage. If for a process stage p the item i must be placed on 
machine k, we require 


n(k) 


(1) v: xt) = 


j=1 


That is, the item i must appear in some ordered position on machine k; the restriction that 
x) be integer-valued implies that the item will be scheduled once, and only once, on machine 
k, If there are Q machines listed for process stage p and the item must be placed on every 
one of the Q machines, then we have a set of Q equations of the form of Eq. (1), one for each 
such machine. A similar statement holds for the relations below. 

If for a given item i and process stage p the item must be placed on one machine out 
of a group of machines ky, Ko, mole ko, we require 


n(k,) n(ko) n(ko) 
k,) (ky) (k 
(2) x," + x" Hoot x2 oo Bg 


j=1 j=1 j=1 


Finally, if for a process stage p the item i must be placed on any of q machines out 
of a group of machines ky, Ko, ... , Kg, we require* 


n(k) 

(k,) 

(2a) >. Xj 1 
j=l 
n(Ko) 


yo 
J 


j=l 


°Q 


6; +59 +--+ +59 = Q- a4, 





*The general form of the constraint suggested is given by Dantzig [7]. 





134 H. M. WAGNER 


where each 5. is restricted to be 0 or 1. Equations (2) and (3) and the integer restriction on 
6 ensure that item i is processed on exactly q out of the Q machines. 

The second set of constraints guarantees that no more than one item be assigned the 
j-th ordered position on a machine. For each machine k, we require 


(4) z; i Tak woe | Om 
ie N(k) 


It should be noted that (1), (2), (3), and (4) and the condition that the x) and 5, be 
integer valued are sufficient to restrict these variables to the values 0 and 1, i.e., no 
additional upper bounds are needed. 

Reviewing the above constraints, we recognize there is yet no guarantee that (1) an 
item be scheduled for process stage p and completed before it is scheduled for process stage 
p +1, or even that (2) an item not be scheduled to be processed on two or more machines 
during the same time. Consequently we need a set of (linear) relations implying that the 
production schedule observes the process-stage restrictions and does not call for the item 
being placed on more than one machine at one time. Toward this end we find it convenient to 
suggest a "shorthand" notation and to give explicit relations for n{X), 

Let 


k k k ‘ 
(5) Ty) * ;> 1(*) ~ o Ss. , eed. 
ie N(k) 


Given (4), T (k) represents processing time of the j-th ordered item on machine k. Then 
for each machine k it may be verified that 


(6a) nf) = s(*) 


(6b) n{*) = Tx) + T x} +eeet Tx) + s(*) + s{*) eu2ce aX) . j = 2, 3, ... , mk) 


Equation (6) states that the item in order-position j on machine k commences processing at a 
time equal to the sum of the manufacturing and idle periods accumulated from the initial com- 
mencement of production. 

First we consider the restriction that for a particular item i any machining taking 
place in process stage p + 1 may commence only after all machining is completed in process 
stage p. We suppose in the i-th row of matrix I machine ky appears in the entry for column 
p and machine Ko appears in the entry for column p +1. Suppose further that 
(ky) (kp) ei be a 
xjr = 1 and Xin = 1, i.e., item i is scheduled in order-position j' on machine ky, and 
in order-position j"’ on machine Ko. For the schedule to be feasible 


(c,) ey) 5) C) 
hy: + tj ij’ s i" 


(7) 
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Given our hypothesis about the specific order of production for item i, (7) guarantees 
that machining on (Ko) is not commenced until machining on (k,) is completed. But we can- 
not add (7) in its present form as a constraint, since it requires too much, viz., that the 
starting time of the j"-positioned item on machine Ko never precedes the finishing time of 
the j'-positioned item on machine ky; we want (7) to hold only whenever item i happens to be 
the item scheduled in both of these ordered positions. Therefore we merely need to require 





" ed Od Oe) Oe) to) 
i" ,-_— = = i ae ill i 
where M is a large positive integer and nf) is evaluated by (6). Observe that (8) is a binding 
constraint only if a? = ov = 1 
J 1) 
In general, for item i, each pair of process stages p and p +1, each ordered couple 
(machine k, in process stage p, machine ky in process stage p + 1), and j'=1,2,... »n(k,), 


<a n(Kp), we have a constraint of the form (8).* Although the total number of such 
constraints for the model is obviously enormous, only a relatively few of the relations will be 
binding in any solution. This fact suggests the technique of treating relations of the type (8) as 
secondary constraints [4, 8, 14]; in other words we might attempt to solve a scheduling problem 
by a series of trial solutions in which constraints are introduced only as needed to eliminate 
infeasibilities. 

Secondly we consider for item i the restriction that any machining taking place within 
process stage p must be on only one machine at atime. We suppose machines ky and Ko 

ae ; (Ky) 9) 
appear in the i-th row and p-th column of matrixI. Under the assumption Ar = Ayer = 1 


? 


we require that one of the relations below must hold in order that the schedule be feasible: 


k 
(9a) a? +t 


(iy) Cy) a) 
j 


i Xj" < a. , 


(5) ) (Ky) 
Din + ti se Sh j 


(k,) 


(9b) h., 


Analogous to the reasoning in (8), we want either (9a) or (9b) to hold only whenever item i 
happens to be the item scheduled in both of these ordered positions. Our linear constraint then 
is the pairt 


k k k k k k 
(10a) ny ) + ‘ ) x) < hyn + m( - a) + m( - i? +6M, 


k k,)  (k k k k 
(10b) hyn? + ‘ 2) x2 . ry + m1 ‘ | +M(1- a | - M(5- 1), 


where 6 is either 0 or 1. 


*As the reader may verify, we need not consider constraints between anything more remote 
than consecutive pairs of process stages. 


+The general form of the 6 constraint suggested is given in [7]. 
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In general, for item i, each process stage p, each couple (machine kK, in process 
stage p, machine kp in process stage p), and j'=1, 2,..., n(k,), _ ae n(Ko), we 
have a constraint of the form (10).* - 

We finally come to the question of the objective form. Let h represent the earliest 
point in time at which processing on all items has been completed according to a given trial 
schedule; we desire to find that schedule which minimizes h*, If, for example, all items must 
be "finished-off" on machine m, then the optimizing form is simply minimize (n™ m) , T xi™) ), 
where (6) is used to evaluate ni™), If a "finishing" machine does not exist and there is no a 
priori reason to believe that a particular machine will be the last one in operation, then we 
must allow for the possibility of any machine being the last in operation. One method of 
solving the problem is to minimize h*, subject to the additional constraints 


(11) ny + Txhy ch Ts Ss oP 


where, as in (8) and (10), the term on the extreme left hand side of (11) is evaluated by (6).+ 

To summarize, we have formulated a scheduling model capable of handling, theoreti- 
cally, a wide class of machine-sequencing problems. In spite of the fact that the total number 
of constaints to be satisfied is staggering, most of them are inoperative, and the technique of 
secondary constraints may prove useful. In addition to the complexity of the problem due to its 
size alone, a further computational difficulty is introduced by the requirement that an integer 
linear-programming algorithm of some sort be employed. Nevertheless, the demonstration 
that there does exist a finite procedure which monotonically searches for an optimal solution‘ 
affords some encouragement that eventually a practicable algorithm may be developed. 


Il. THE CLASSICAL PROBLEM OF THE m-STAGE PROCESS 

A well-known special case of the above model is the problem: Given n items, each to 
be processed on m machines, the order of processing being specified as machine 1, machine 
2,..., machine m, find the sequencing of items on the machines that minimizes the time at 
which the last item is completed on machine m [2, 3, 10, 12]. 


Now all items must go through m process stages, machine k is specified for process 
stage k, and n(k) = n. Constraints (1) and (4) for each machine k = 1, 2, ... , m become 


*As the reader may verify, we need consider only all pairs (k,,k z) to ensure feasibility 
within process stage. 


+tAn alternative device would be to define a ficticious "finishing" machine with process time 
equaling zero. 

tA comment is called for in explaining the meaning of a "monotonic" search. The modus 
operandi of most extant integer programming algorithms is to ignore the integer constraints 
and to find an optimal solution to the linear-programming model by a standard technique; ifa 
nonintegral solution is thereby obtained, constraints are added one by one whicheventually force 
the enlarged problem to yield an optimal integer solution to the original model. Consequently 
the value of the optimizing form is first (monotonically) brought, say, to a minimum, and then 
monotonically increased by the addition of new constraints. Similarly in our model, if both the 
integer restrictions as well as the feasibility constraints are added as needed, the optimizing 
form would first be minimized and then increased monotonically as feasibility is established. 
The important claim for a finite monotonic procedure, of course, is that in most problems it 
eliminates having to enumerate completely all possible solutions in orderto findan optimal one. 
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n 
2. 


j=1 
i ar 


Notice that for each value of k, the relations (1') and (4') describe a standard “assignment- 
problem" model. 


By the use of (6) and (1') evaluated for machine m, the minimizing form reduces as 
follows 


(12a) ni) + Tx'™) - Txi™) + Tx") 4 + Tx™ 4 5 (m) + sim) +eeet s(m) 


im) (xin) + xim) Dsante xm) ee ¢(™) (xm) + xim) +eeet x(m)) + 


(12b) + s\m) teeet s(™) 


(12c) = constant + s\m) + s(m) wet s(m) . 


In (12) we have the commonly recognized relation [2, 10] that the optimizing problem may be 
stated as minimizing the total amount of idle time on machine m. 

We may, of course, use constraints of the form (8) to ensure a feasible processing 
sequence. But we may also formulate another set of restrictions, based on this model's 
special assumptions, which encompasses many (although not all) of the restrictions implied by 
the full set of relations (8). To do this we define nonnegative integer valued variables 


ult) = time elapsing between the completion of processing the 
item scheduled in order-position j on machine k and 
the start of processing the item scheduled in order- 
position j on machine k +1. 


We know the item scheduled in order-position j on machine k + 1 cannot commence 
processing before the item scheduled in order-position j on machine k is completed; for, 
otherwise, there would be a violation of the restriction that the machine order of production 
is the same for every item.’ Note that the previous statement does not imply that the item in 
order-position j on machine k is the same item as that in order-position j on machine k + 1; 
the restriction is solely one of the time interdependence among process stages. 

Explicitly for period j(>1) and machine k(>1) we have 


(k) _ ,(k-1) (k-1) | ul 1) 
(13) hj = hj + TX 


506275 O-59-7 





138 H. M. WAGNER 


That is, the time at which the processing of the item at order-position j on machine k begins 
must not be earlier than the time at which the item at order-position j on machine k - 1 is 
completed. 

We may also write (6) in the form 


' k k k k 
6" 10 = nls al) + of). 


Subtracting (13) from (6'), we get 


—- pik) _ _(k-1) (k) _ (k-1) (k) _ ,(k-1) 
0 = a hy + TX; 1 TX + sj uj ‘ 
Subtracting (14) evaluated at j' and k' from (14) evaluated at j'+1 and k' and using 
(6'), we derive the constraints 


(15) 0 = Tat) . pg), gD te'-1)  e'-1) =) 
= j' j'+1 j' 7 j' j'+1 . 


er sem , 8 * fei. 


It may be verified that the (m-1)(n-1) relations in (15) comprise a modified transporta- 
tion problem; each variable appears in no more than two of the relations, although the — 
have coefficients other than 1 or -1. 

The structure of the model [2,10] also permits us to set 


(16) a1) = uf) = uD) = 0 


without loss of optimality. 

Clearly (15) imposes restrictions only on the timing of processing the item scheduled 
in order-position j, but does not refer to any timing constraints on the particular items them- 
selves. Consequently if we were to use a secondary constraint technique to solve a specific 
problem, we would need certain of the relations (8) in addition to those in (15) to reach a fea- 
sible solution. 


IV. THE CASE OF THREE MACHINES 

A particularly interesting case of the model in Section III is with m= 3. This class of 
problems is distinguished by the fact that, without loss of optimality, we may confine our 
attention to schedules which sequence the n items in the same order on all three machines 
[2,10]. Mathematically, this means that we may restrict the values of the unknowns by the 
relations 


et 28. coe OO - he 2 Since Rs 


As a consequence, we need only impose (1') and (4') for k = 1, comprising a system of 
2n relations and n“ unknowns (in the format of an assignment problem, as we noted above). 
Similarly, (15) may be rewritten in terms of =, yielding a set 2(n-1) equations. 
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It may be verified that the minimizing form can be written as 


minimize T x() + 7 x() + s(3) + s{3) +eooet s(3), 

Finally we note that, in this special case, the only timing restrictions necessary are 
those imposed by (15), i.e., there is no need for any additional constraints of the form (8). 
This conclusion follows from the observation that now there is no ambiguity between "'the item 
scheduled in order-position j on the machine k" and the j-th item in the scheduling sequence. 

Thus we have derived a fundamental system of the order of 4n equations which proba- 
bly can be solved* for n < 25 on high-speed computing machinery; but it of course remains 
questionable whether or not such a computational proposal for finding an optimal solution is 
economically sound. 
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Renewal theory is used to derive the stationary distribution of 
inventory present at a stocking point characterized by stationary 
compound Poisson demand, independent and generally distributed lead 
times, impatient customers, and the use of a two-bin ordering rule. 
Inventory-associated costs are evaluated, and a numerical example is 
given. 











INTRODUCTION 

Some variation of the ''two-bin"' system of inventory control is frequently used in 
warehouses to control finished goods stocks and in factories to keep stocks of raw materials 
and partially finished goods at suitable levels. This system usually operates about as follows: 
The total stock of an item is divided (possibly physically, possibly only conceptually) into two 
parts: the "reserve stock" and the remainder, which can be called the "ordinary-usage stock." 
The reserve stock can be thought of as occupying a "second bin" which is kept locked until 
customer demands have entirely depleted the ordinary-usage stock. When that customer 
demand arrives which actually exhausts the ordinary stock, an order for more stock is placed 
with a supplier, and the contents of the "second bin" are made available for issue. Actually, 
the customer demand that eventually depletes the ordinary stock will almost certainly also 
require a dip into the reserve stock, so not quite the entire reserve stock is available after 
the order goes in. The purpose of the reserve stock is, of course, to satisfy customer 
demands during the "lead time,"' where this is the time that elapses between the placing of the 
replenishment order and the arrival of this order. In practice, no "second bin" need exist 
physically, but when the total stock on hand drops below a preassigned level, called the "order 
point,’ an order is placed for replenishment stock. This replenishment, which will henceforth 
be called simply an "order," arrives one lead time later. 

In order to use such a system, the controller must decide on the amount of reserve 
stock to be kept on hand (the order point) and the quantity of replenishment stock to be ordered 
when an order point is cut. It is reasonable to expect that management responsible for inven- 
tory control will want to maximize returns from inventory, whether this aim actually is made 
explicit or not. Under static conditions of price and production technology, maximizing returns 
is equivalent to minimizing the costs associated with carrying a certain general level of inven- 
tory. Assuming that conditions are static enough to make cost minimization a reasonable 





*Manuscript received April 22, 1957. 
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objective, an order point and order size can be determined so as to attain this objective when 
a particular type of demand pattern is anticipated. How the minimizing values or order point 
and order size can be found is discussed in this paper whenever a particular type of random 
demand is postulated. The optimizing order point and order size depend, of course, as much 
on the demand pattern as on the costs. Specification of both these inputs will be an important 
part of any application of the theory, but one which is not considered here. 

It is frequently reasonable to assume that there are three costs associated with carry- 
ing inventory stock: (a) carrying costs, (b) ordering or restocking costs, and (c) stockout 
costs. These somewhat self-explanatory terms are defined in Section 5. The expected 
(average) yearly value of each of these costs is affected by choice of the order point and order 
size, and their sum can be minimized by such a choice. When replenishment stock is available 
instantly (lead time is zero) no protection stock is necessary, and a simple formula can be 
derived for the order size that minimizes the sum of costs (a) and (b). This "economical-order 
quantity" is well known [1]. The effect of the nonzero lead time is to make necessary a certain 
amount of reserve stock to keep stockout costs sufficiently low. That the optimum levels of 
protection stock and order size cannot be determined independently is also known, but the 
relationship seems to have been explored very little. In this paper a formulation of the prob- 
lem is given when nonzero lead times exist, and although optimal levels of order point and 


order size do not emerge as simple formulas, a moderate amount of numerical work should 
yield them. 


Considerable important work has been done previously to develop optimal strategies 
for reordering stock when certain types of costs are assumed [7-12]. The objective of this 


paper is the more modest one of describing inventory fluctuations when a particular type of 
reordering strategy is in effect: the so-called "two-bin" rule. Finally, some possibly 
interesting system costs are discussed and evaluated in terms of the earlier results, making it 
possible to choose the parameters of the "two-bin" rule (order point and order size) optimally. 
A numerical example is given. 

The methods of this paper are those useful for analyzing many types of simple stochastic 
processes encountered in applications, such as in the theory of queues and waiting times [2], of 
renewal theory [3], and of Geiger counters [4]. Renewal theory also has been applied to inven- 
tory studies [7,13]. Thus, even though the specific model analyzed in this paper may not be 
applicable in a particular situation, it is hoped that the same methods may perhaps be used by 
the reader to construct one that is. 


1. OUTLINE 

In this paper is derived the stationary distribution of inventory on hand when a "'two- 
bin'"' system of inventory control is in use. The stationary distribution is then used to evaluate 
an inventory cost function, which in turn can be minimized by choice of an order point and an 
order size. Demand is described by a random process that is itself stationary. This means 
that the probability distribution of the number of items demanded in a time interval of fixed 
length does not change if the initial point of the time interval is shifted about on the time axis. 
Another way of putting this is that demands during nonoverlapping periods of time, such as 
successive months, are independent and identically distributed random variables. 

It is further assumed that demands that cannot be satisfied immediately with stock on 
hand do not wait (there is no back ordering), but that their disappearance does not mean a drop 
in the expected number of demands per unit time, i.e., there is no loss of "good will." This 
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assumption represents a situation in which there are effectively many customers with short 
memories whose average demands for the item being stored change little with time. 

The effect of a nonzero lead, or delivery, time (the time between crossing an order 
point and the arrival of the stock orders) is included, but the ordering procedure analyzed 
permits only one order to be out at atime. This restriction is probably not practical in case 
the lead time is long enough, so that many demands are likely to occur in its duration. 

This model applies strictly only under the mentioned restrictions, but practically it 
may be useful even when, for example, occasional extremely large demands are received that 
management feels should be filled, even though doing so means temporarily abandoning the 
system. The existence of these orders might be recognized by simply suspending ordinary 
operations when they arrive and concentrating upon them alone or, depending upon their 
frequency, it might be best to simply mix them in with the rest of the demands. The problem 
of when executive decision should override system decision is an important problem that has 
not been treated here. 


2. ASSUMPTIONS 
Specifically, we assume that demand is a time-independent compound Poisson process 


with rate parameter \, so that the probability of receiving exactly n demands for the item in 
question in a fixed time t is 


yt (at)™ 
(2.1) ert mae 


The average number of demands is ) per unit time, thus ) is the demand rate. In order to 
estimate \ in practice it is simply necessary to count all demands for the item in question 
occurring in a unit time interval. If the time interval chosen includes, or is included in, a 
period in which the supplier is out of stock, the count must include the number of unfulfilled 
demands as well as the number of those for which stock is available in order to be a valid 
estimate of the demand rate. 

It is further assumed that demand sizes are independent random variables, with the 
common distribution function 


(2.2) 


Prob. {z < z} 


and E (2) 


where Z is the size of a demand chosen at random. Exponential distribution (2.2) was selected 
to represent demand size variations partly because empirical study has shown it to be not 
unreasonable in some cases and partly because of its desirable mathematical properties. That 
any real number is a candidate for a demand size is a convenient and noncrucial assumption. 

If stock were initially at some level Ij, then inventory would tend to decrease more or 
less steadily according to some step function ("sample function") such as that shown in Fig. 1. 
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The times between successive downward jumps 
are independent random variables with the dis- 
tribution function 








Figure 1 This follows from (2.1). The heights of the jumps 
are independently distributed according to (2.2). 
The assumption of a process in "continuous time" is again merely a noncrucial convenience. 
In order to satisfy demand on a continuing basis it is necessary to replenish stocks 


periodically. Therefore, we introduce the following replenishment rule (two-bin or order- 
point method): 





Whenever the amount of inventory of a type of unit on hand, I, falls below the 


fixed level, r, the order point, order an amount, R, the order size, from a 
supplier. 


The time between the placing of an order and its receipt, the lead time, is generally not 
equal to zero and may exhibit random fluctuations. We assume that successive lead times are 
independent random variables with common distribution function U(z), so that if L is a lead 
time chosen at random, 


Prob. {i <2} = Us). 


U(z) will be assumed to be absolutely continuous, so that 


Zz 
U(z) = l u(x) dx, 


where u(x) is the probability density function of lead time. We also assume that the mean 
(average) lead time 


@ 
eq) - [ z dU(z) 


is finite. Neither of these assumptions restricts the applicability of the theory. Lead times 
are assumed to be not only mutually independent, but also independent of the state of the 
inventory, so there is, for example, no expediting. The entirety of an order arrives at the end 
of the lead time in the present model, although in practice large orders might reasonably be 
expected to be filled in installments. 

From the foregoing remarks, it is clear that for fixed levels of r and R(R>r) there 
will be a certain probability that the amount of inventory on hand at time t does not exceed the 


level i, providing the inventory was equal to iy initially. Thus 


Fo (i, t) = Prob. {h silly = io} ; 
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Moreover, the stationary distribution‘ 


F(i) = lim Fo(i, t) 
t 


—>o 


can be shown to exist and to be independent of the initial conditions ip: Both Fo (i, t) and F(i) 
are functions of the demand parameters \ and w, the lead time distribution U(z), and the 
replenishment parameters r and R. 

Having the stationary distribution, three important measures of the effectiveness of the 
inventory system can be calculated in terms of the above parameters: (i) The stationary 
probability that there is no inventory present (of being "'stock-out") F(0). (It will be shown 
that F(i) is interpretable as the expected fraction of a long period of time during which the 
inventory does not exceed the level i, so F(0) is the expected fraction of a two or three year 
period, say, during which the supplier is "stock-out".) (ii) The average level of the inventory 
on hand 


~0 
E(I) = if i dF(i). 


Finally (iii) the average time between the placement of successive orders for more stock. A 
discussion of each of these measures and their associated costs will be given in Section 5. 

Use of the replenishment rule described above is not a unique or necessarily optimum 
way of controlling inventory. Whether or not following a particular replenishment rule yields 
truly optimum (minimum) total inventory cost over a period of time depends upon the costs 
attached to ordering, carrying inventory, and being out of stock, as well as upon the fluctuations 
of demands and lead times. Many writers [8-12] have described attacks on the problem of 
determining truly optimal replenishment rules, given demand and lead time statistics and the 
cost structure. It turns out that in many cases of practical interest optimal replenishment 
rules are exceedingly complicated. In practice one is frequently forced to select a subclass of 
the class of all possible rules (e.g., the subclass of two-bin rules) and to optimize within that 
subclass (e.g., by choice of R and r). The latter approach is taken in this paper. A numerical 
example appears in Section 6, in which the results obtained here are compared with those 
obtained with the more classical methods [1]. 


3. THE STRUCTURE OF THE STATIONARY DISTRIBUTION OF INVENTORY 

When the order point is cut, i.e., stock passes from a level above the order point to 
one below it, the order cutting it reduces the stock to some level below r. Let the random 
variable Y denote the magnitude of the "overshoot." A property of the distribution (2.2) that 
proves to be of great convenience is that the overshoot Y has the same distribution that the 
size of a randomly selected demand has: 


(3.1) Prob \y < y| = 





tF (i) is not the stationary distribution of a Markov process. Our inventory process is in 
general non-Markovian. 
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Furthermore, Y can be shown to be independent of the time when the order point was last cut. 
This implies that the times between successive cuts of the order point are identically dis- 
tributed, independent, random variables. In fact, given that the order point is cut at time t 
with an overshoot Y, no further information about events occurring prior to t has any bearing 
on events which will occur after t. These facts, which are true in general only for the expo- 
nential distribution of demand sizes, greatly simplify subsequent calculations and make it 
possible to determine stationary distributions which are informative about the problem of 
inventory-control policy. 

In order to obtain the stationary distribution of inventory present, F(i), it is convenient 
to find first the distribution, F(i, t), of inventory present at time t, given that initially (i.e., at 
t = 0) the inventory level had just dropped below r, the order point. Let N(t) be a random 
variable denoting the number of times the order point is cut in the time interval (0, t), 
excluding the intial cut. Then clearly 


F(i, t) = Prob {I < i} 
(3.2) 
= G, (i, t) + Go (i, 9%. 


where 
(3.3) G, (i, t) = Prob {¥ i, N(t) 


(3.4) Go (i, t) = Prob {hy i, N(t) > 
In order to find Gy (i, t), we write 


t 
(3.5) Gy (i, t) = { G, (i, t-7) dQ(r), 


where 


dQ(rT) Prob | N(r+ 47) > N(r)} : 


the probability that an order point is cut at least once in the interval (7,7 +d7T), Q(T) being 
Prob { N(T) > 1}. As will be clear later, there is actually a "reorder density" q(T), so 


dQ(r) = q(7)d7 = Prob | N(r+d7) ° N(7) +1} 


as dr>0. 


Now let X be the reorder time, the time between two successive cuts of the order point, 
and let 


(3.6) Q, (x) = Prob {x < x} ‘ 
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This distribution is absolutely continuous, with density qy (x), and a finite first moment, the 
expected (average) time between successive cuts of the order point, given by 


foo) 


E(X) = | x dQ, (x). 


Now, if initially the order point has just been cut, 


© 
(3.7) Q(x) = >; qn" (x) , 


n=1 
where the asterisk indicates the convolution operation: 
1* 
( 


Q™ (x) = Q,() 


x 
ate = | oP - 2) aQe). 


Formula (3.7) simply expresses the fact that an order point is cut at some time before the 
time x if this happens exactly once, exactly twice, ... The sum of the probabilities of these 
mutually exclusive events gives Q(x), since if X, is the time between the (n- 1)t and nth 
cuts of the order point, the sum KX, +Xo+--- + Xx, is the time to the nth cut, and 


n* 
Prob {Ky + Xp tee +X < x = Q; (x) . 


Combining equations (3.2) to (3.5), we have 


t 
(3.8) F(i, t) = G, (i, t) + | G, (i, t -7) dQ(T). 


It is shown in Appendix I that the limit as t tends to infinity of F(i, t) exists, is 
independent of the initial conditions, and is equal to 


© 
i} G, (i, t) dt 


E (X) 





(3.9) F(i) = 


This is the form of the stationary distribution we are after and which will be constructed in the 
next section.t 





+The formula (3.8), with what follows, provides the distribution of inventory at any finite 
time in future. The latter can be used to evaluate discounted costs and to optimize ''dynamic" 
policies. This will not be discussed further in this paper. 
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4. CONSTRUCTION OF THE STATIONARY DISTRIBUTION 

The distribution functions Q, (x) and G(i, t) that are the necessary ingredients for the 
stationary distribution (3.9) are both constructed from the basic demand distributions (2.1) and 
(2.2) by making use of the reorder rule and the hypothesized behaviour of the customers. 
Proceeding towards these goals, let Vv; be the total quantity of stock demanded by customers 
in a period of length t. If t is measured from the instant the order point, r, is cut and the 
overshoot, Y, is included, then Vt is a random variable with distribution function given by 


© (at)? v/w 
(4.1) H, (t, v) = Prob {V, < v| = ze: gens ee | 


n! 
n=0 


The integral 


v/w n 


x * 
eX — ax = wirl) (v) = Prob [Y¥+Vp+Vy+...4V cv}, 
0 n! n 


Vj being the size of the jt demand. Multiplying wntl)* (v) by the probability that n demands 
take place during the interval of length t gives the desired probability distribution. 

Note that, upon subtraction from unity, Hy (t, v) is also a distribution function, now with 
respect to t, v being considered as fixed. This is to be interpreted as follows: Let Ty be the 
time after the initial cut of the order point at which total demands first exceed the fixed level 
Vv; again the overshoot is included in v. In order for ¥e >t, it is necessary and sufficient that 


Vt <v, so 
(4.2) H, (t, v) = Prob {V, < v = Prob {ry > t} . 
and 


(4.3) 1 - H, (t, v) = Prob {Ty < t| . 


Considered as a distribution function with respect to v, t being fixed, H, (t, v) is absolutely 
continuous (possesses a density) for all v, 


= (at)? 
(4.4) h(t, v) = >: gt e 


n! 


n=0 


As a distribution function in t, 1 - Hy (t, v) possesses a density for t > 0, but at t = 0 there 


is a finite jump of height e“¥/ ®. This corresponds to the fact that the demand cutting the 


order point may overshoot by an amount v or more, and this happens with probability o*/ a! 
It also will be necessary to have the distribution of the amount of stock demanded in 


time t, Vi» when there is no initial overshoot. The same reasoning that was outlined above 
shows that, in this case, 
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v/w -1 


at, Ft grat HM : 

- -xX 

(4.5) Ho (t, v) = Prob {v, < v} e + 2, e 7 l e (n-1)! dx. 
n= 


Again we have 
(4.6) Ho (t, v) = Prob {vy < v} = Prob {r, > t} ‘ 


Considered as a distribution function in v, Ho (t, v) has a jump of height s* t at v= 0, cor- 


responding to the probability that no demand occurs before time t, but for v > 0 it possesses 
a density. The distribution function in t, for fixed v, has a density in t for all t >0. 

We now construct the distribution function Q, (x) and its mean, E(X). Observe that X 
is the sum of the lead time, L, and the passage time, T. The latter is the time that elapses 
while the inventory reduces from its level just after an order is received to the order point. 
Note that there is a dependence between T and L, for if by chance L happens to be long, more 
stock will tend to be demanded than would be if it were short; consequently the level of inven- 
tory would tend to be low just after the new order arrives, reducing the time to cut the order 
point again. 

If v. = v is the stock used during a lead time of length z, then the stock level just 
after the order arrives is R+r - v, and the stock that must be used after the order is 
received to cut the next order point is (R+r- v) - r=R-v. Thus, given L = z, the proba- 
bility distribution of the passage time, T, is 


ial P(t, z) = Prob {r, < t| 
r 


= | {1- Hy(t, R-»)} d, H,(z, v) + {1-Hy(t,R-»)} {1 - Hy (2, )} ° 


The integral is the sum of the probabilities, respectively, that the passage time is less than or 
equal to t, given that the stock to be used up just after the order arrives is R - v, multiplied 
by the probability that the stock to be used up is indeed R-v, given the length of the lead time. 
This part of the sum represents the situation when no stockout occurs. The last term repre- 
sents the probability distribution of the time to cut the order point again, given stockout 

(Vv, > r), times the probability stockout occurs in the leadtime. Finally, by multiplying by the 
probability that L = z and adding over all z, we obtain 


[> 6) 
(4.8) P(t) Prob {7 < t} = [ P(t, z) dU(z). 


Since X= T +L, 


x 
(4.9) = Prob {x < x | = if P(x - t) dU(t). 
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Since the stationary distribution F(i) depends only upon E(X), we shall evaluate this mean as 


explicitly as possible. In Appendix II the Laplace transform of the distribution of the passage 
time, T,, is derived, and from it, by differentiation, the mean passage time. The result is 


- 1 Vv 
(4.10) E (TY) = | td H(t, v) *. * So 


Using this with (4.7), we have 


E(T!z) = 


io) 


- eye oe -r»z Lz)? 
“5° - 


AW Aw n! 
n=0 


By multiplying this by the probability that the lead time takes on the value z and integrating, 
we obtain 


Since the expectation of the sum of two random variables is the sum of the respective expecta- 
tions, whether or not the variables are independent, 


(4.13) E (X) E(L) + E(T), 


where, of course, 


re) 
E(L) = l z dU(z). 


Construction of the distribution function G(i,t) naturally divides itself as follows: 

Case (i): O<i<r 

In order that I <i, stock must have declined by an amount at least equal to r-i in 
time t, the time since the order point was last cut, and the lead time must not have terminated. 
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Thus 
(4.14) - {1-H tt, r-a} {1 - ve}. 
For example, 
(4.15) 
and 
G0, t) = {1- H(t, 9} {1 - vw}. 


Case (ii): r<i<R+r 

It is convenient in this case to compute the probability that the inventory exceeds i. 
This happens if total demand after the order arrives does not exceed R+r-i-v, v being 
the demand during the lead time. Thus 


G(R +r, t) - G(i, t) = 1 - G(i, t) 
(4.16) 


tos 


t [ 2 
= I dU (z) f falta, Res-to¥) A Hye, g(x, Ro) (1-H) | 


By referring to (3.9), it is easy to show that 


we 
(4.17) if G(R+r, t)dt = E(X). 


This indicates that F(i) obeys one of the rules for well-behaved distribution functions, namely 
that F(o) = F(R+r) = 1. Substantiation follows from (4.16) by setting i = r, changing the 
order of integration, i.e., integrating first on t, integrating by parts once and referring to 
Appendix II, and finally integrating on z. Without going through the details, this verifies that 


io) 


(4.18) | {G(Rsr, t) - G(r, »| dt = E(T). 


Integrating (4.14) by parts once verifies that 


co 


(4.19) i. G(r, t) dt = { {1- vw} dt = E(L). 


Addition completes the result. The reversal of the order of integration in the first step above 
is justifiable by Fubini's theorem, or, more simply, by reference to Titchmarsh [5, p. 54]. 
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By reversing the order of integration and integrating once by parts—the same manipu- 
lations that lead to equation (4.17)—we have 


1 R+r-i-v 


co R+r-i 
(2 Bere) dH, (v) R<i<Ré+r 


| G(R+r,T) - {| Gti,7)| dr if 


4 a. R+r-i-v 
heme mes dH, (v) 


.( + “ (1 - Hy (x)) , 


where for convenience we have written 


H, (v) 7 | H, (t, v) dU(t), 


the (unconditional) distribution of the stock demanded during the lead time. 
Finally, from the use of (3.9), the stationary distribution of inventory stock on hand is 


R+r-i R ; 
nm - [ BEY) oo) 


r 





ace —, R<i<R+r 
E (X) 


r . . Ps 
E(X) - { (2. Betty) d,, H, (v) (2 Bt) a ‘ i,t 


a ; r<i<R 
E (X) 





@ 
E(L) - | H(t r-9{1-veo| dt 
. E (X) 





Equations (4.21) suggest the following interpretation: F(i) is the average fraction of a 
reorder time during which the inventory level is below i. This is easily seen by observing 
that the integrand of the first two expressions 








is just the expected time taken by the inventory to pass from the level R+r-v immediately 
after an order arrives to the level i(i < R+r-v). Thus, given v, the amount of inventory 
removed by demands during the lead time, E (T; ,) is the expected time spent by the inventory 
above the level i during a passage, and hence a het time. Multiplying by the distribution of v 
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and integrating gives the unconditional expected time spent above i, and subtraction from the 
expected length of the lead time gives the expected time spent below i. Similarly, rearranging 
the integral in the numerator of the last expression gives 


io.8) fe) [8 
(4.22) | G(i,r) dt = | {1 - u(r)} és | Hy (7, r-i) {1-u(r)} ar. 


The first integral is the expected length of a lead time, from (4.19). The second is the 

expected time to use an amount of stock equal to r-i when the lead time has not elapsed; for 

H, (t, r-i) j1 - U(t)} is, from (4.2), the probability that the time to use an amount of stock 
equal to r-i exceeds t and that the lead time also exceeds t, by the independence assumptions. 
An integration by parts establishes the result. The difference between the two integrals is the 
average time the inventory spends below the level i(i < r) during a lead time, or, equivalently, 
during a reorder time. It follows finally that the stationary probability F(i) gives the expected 
fraction of a "long" time for which the inventory level does not exceed i. This comes, for 
example, from (A-3) of Appendix I. 

Differentiating equations (4.21) gives the probability density function of the stationary 
distribution, where it exists. There is a nonzero probability that the inventory will take on the 
value zero, and hence a nonzero probability that it will take on the value R just after an order 
arrives following a lead time during which stockout has occurred. Representing the stationary 
distribution in terms of densities, we have 


=. H, (R+r-i) + + hy (Rer-i) 
r 


_ AW 
E (X) 





1 eet 
— (1 - Hy (x) 


E (X) 





1 
Ne 
E(X) ’ 


4 hy (t, r-i) {1 - vie} dt 
E (X) 





i. H, (t, r) {1 - vie} dt 


E (X) 





’ 


where we write 
d F (i) 
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if this probability is not zero. Notice that for part of its range (r< i < R) the stationary 
distribution is uniform, a result that is not surprising when the inventory process under 
discussion is compared with the random walk over a finite interval. 

All moments, and in particular the first, of the stationary distribution follow by 
integrating (4.23) over the range. Thus 


(4.24) = i* dF(i). 


The resulting formulas are not simple, in general, so they will not be written here. The 
expressions are easy to write, but not necessarily to evaluate, in any interesting special case. 


5. INVENTORY COSTS 

The stationary distribution permits us to evaluate a particular inventory control 
system, i.e., one characterized by known parameters \ and w, known lead-time distribution 
U(z), and fixed R and r, from the standpoint of costs. Since the values of R and r, in an 
actual situation, are subject to the control of the supplier, they can then be selected to minimize 
the total cost of carrying inventory when the type of two-bin control assumed here is in use and 
customers do not wait. 

Costs will be assumed to come from three sources: (a) carrying charges, such as 
warehouse rent, insurance, etc., (b) ordering costs, such as machine set-ups, order processing, 
etc., and (c) stockout costs, in this case the loss of the immediate business of customers who 
happen to find the warehouse empty. Our stationary distribution allows us to evaluate each of 
these separately on a long-term basis by arguing as follows: let t be a "long" period during 
which customers are supplied from inventory in accordance with our assumptions. Then F(i) 
is the average fraction of the time t during which the inventory does not exceed level i, so 
td F(i) is the time during which inventory is "at" level i. The carrying charge associated with 
the inventory at the level i is proportional to the level i and also to the time that inventory is 
actually at that level. Summing over all levels, the carrying cost in dollars for atime t to 
terms of order t, is 








R+r 
t ¢, { i dF(i) 
0 


t m, (R, r), 


fe 
where m, (R, r) is simply the mean of the stationary distribution described in the last section. 
Other factors held fixed, the value of m, (R, r) will increase as R and r increase, so carrying 
costs will increase with increases in the order point and order quantity. Parenthetically it 
should be remarked here that the carrying charges evaluated by (5.1) are not the only ones that 
can be considered. For example, it might be reasonable to allot a charge that varies in accord- 
ance with the maximum possible level of inventory, R+r, for if R+r exceeds present storage 
capacity it may become necessary to build new warehouses or store goods temporarily (at high 
rates) in freight cars. Since the stationary distribution does not depend on costs, any (integrable) 
cost function can be averaged with respect to it. 
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The average number of orders made during the long time t is = according to a 


basic theorem in renewal theory [3]. Allotting a constant cost, c,, to each order, the total 
number of dollars spent on ordering during time t is 


tc, 
(5.2) C. = E(X) 


Since E(X) increases when either R or r are increased, the ordering costs decreases 
with increase in order point or order quantity. The ordering costs considered here are not the 
only possible ones. For example, ordering costs might vary with order size. 

The average time spent in stockout is the average fraction of the time t during which 
the inventory level is zero. This is t F(0). The average number of different occasions of 
stockout during time t is 


(1) t ie 8) wD (x z)" 
— que “AZ — ” 
(5.3) t F’*’ (0) = “ om | e 7 dU (z) { . 


At each occasion of stockout we assume that the part of the customer's demand that cannot be 
filled is lost, as are all customer demands that are made during the remaining lead time. Thus 
the expected total value of the business lost during t due to stockouts is, in dollars, and to 
terms of order t, 


(5.4) C., = tldF Ox, + FY @c., 


where Con is the expected stockout cost per customer. It is sometimes reasonable to set 


-p, 


where w is the average order size and p is the profit to be made from the sale of a single 
item. Like C., the stockout cost decreases with increases of either R or r or both. 

Each of the costs mentioned is the mathematical expectation or mean value of a certain 
random variable. The random variables in question are not independent; but since we are 
interested in their sum, the individual expectations can simply be added, and the total cost 
over the long period t is 


c=2 C€ ++ 


c so’ 


and cost per year, say, is 


' Cc 
(5.6) ¢hea «= . + {ea dF(0) + c,, FO) oh. 
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Formula (5.6) may be used in several ways. If it is felt that all of the assumptions 
outlined earlier are reasonably well satisfied, C' can be minimized by proper choice of R and 
r. The minimizing R and r then constitute an optimal ordering policy for two-bin control 
when customer behaviour is as we have described it. 

Anothev use for (5.6)—in particular the last term—is to evaluate the implied costs 
associated with present choices of R and r. For example, the value of \ Cue d F (0) +C., F 1)() 
for existing levels of R and r estimates the yearly dollar volume of business lost through 
stockouts, or, alternatively, \dF(0) estimates the average number of customers who encounter 
a stockout condition per year if R and r are kept at the present level. Thisfigure may also be 
of use in explaining or predicting fall-off of customer good will. 

Another approach to the inventory policy problem is to fix the average number of stock- 
outs per year, \dF(0)+F 1 (0), by executive decision, thus establishing a relationship between 
R and r. Finally, the sum of expected carrying costs and ordering costs can be minimized, 
subject to the restriction that the average number of stockouts per year be fixed. This approach 
may be a reasonable one, since it allows executive judgment as to the importance of avoiding 
stockouts to influence the decision rule. 


6. DISCUSSION OF AN EXAMPLE 

Investigation of methods for making optimal choices of order size, R, and order point, 
r, are currently underway and will not be described here. Values of R and r that minimize 
total cost rate (5.6) are complicated functions of the parameters of the demand process and the 
costs; in practice, these values can probably best be found by use of a digital computer and a 
combination of analytical and "direct search" methods. 

In this section we shall show that the labor necessary in using our model to find optimal 
values of R and r may sometimes be justified. In order to do this we will: 


(1) Set up a sample inventory control situation, i.e., we state demand parameters, 
costs, and the lead-time distribution function and specify that the two-bin system is 
to be used. 


(2) Obtain order size, R, and order point, r, using a method of Whitin [1]. We term 
this method "classical"; the corresponding values of R and r obtained by its use will 
be designated by subscripts: Ry and Ty: 


(3) Evaluate the total inventory cost resulting from using Ry and Ty: The cost 
function used will be derived from the model developed in this paper. 


(4) Show that there exists a value Ro + R, anda value ro +r,, for which the total 
inventory cost is lower than that obtained by using R, and Ty: 


We shall first present the example, and then discuss it. 
For simplicity we take the lead-time distribution U(z) to be exponential, 


(6.1) U(z) = 1-e 42, 


Making use of this distribution, we evaluate the total cost rate c(R, r). We shall not go 
through the algebra. The result is: 
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2 2 f 

‘ sais eainiadeaa a 
—. 2 a 
1 Ey 


(6.2) rt c(R,r) = 





and 


(6.4) the 


It is possible to make the "set-up," or ordering, cost Cc. a general function of the order size 
R, but this will not be done here. 

In order to obtain simple formulas for the (approximately) optimal values of R and r, 
neglect the contribution made by the lead time to the reorder time, i.e., put Ey = R in (6.2). 
For the moment neglect also the second term in (6.2), and differentiate the remaining expres- 
sion with respect to R. We obtain 


Cs 


C(R, r) = - — + 


1 
r 


=. 
dR 
and by equating to zero and solving we obtain 
(6.5) 

In terms of numbers of units, this is 
(w Ry) ” 


which is recognizable as the classical "economical ordering quantity" formula [1]. 
Having obtained Ry; we now proceed to find Ty: Put E, = R, in (6.2), differentiate 


with respect to r (neglecting the term at! - ihe occurring in the numerator of the right- 


most term of (6.2)), and equate the derivative to zero. The result is 


Ac 
(6.7a) ry, = 3 log, = ; 
a Ri Cc. 
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if determined from this formula, Ty must be taken to be zero in case the logarithm is 
negative. The above derivation is in the spirit of the one outlined by Whitin [1, pp. 56-62] 
for the case of Poissonian demand during lead time. 


An improved approximation to r, can be obtained by including the term = a-e* » 


dropped above. The result is 





1 
(6.7b) rT = ~ log, 


It can be verified that if the distribution function of demand during lead time is Hy (x), then, 
corresponding to (6.7a), we must solve for ry the equation 


a Ry Ce 
(6.8a) 1 - H, (r) = 


A Coy 


and, corresponding to (6.7b), the equation 


Ry co 





(6.8b) S « H, (r) . ; 
ACe + R, co 


Numerical work indicates that the values of R,; and ry given by (6.5) and (6.7a) 
frequently produce total average costs that are close, percentage-wise, to the respective 
minimum costs. It is not yet known how sizeable the differences may be. The situation is 


illustrated by the following numerical example: 


= 4x 10°, 50 (dollars) 


40 (dollars per unit per year) 


Formulas (6.5) and (6.7a) give the "classical" answers 
(6.9) R,; = 200, r, = 183 


in units of the average demand size, w = 5. The cost rate (6.2) can be evaluated at this point 
to give 


(6.10) C (200, 183) = 74,115 (dollars per year) 
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Search reveals that 
(6.11) 
gives a cost rate 
(6.12) C (290, 140) = 71,135 (dollars per year), 


an improvement of slightly over 4 percent. Improvements of such magnitude, if real, may 
sometimes be worth pursuing. Note that the policy of (6.11) differs from that of (6.9) in having 
a larger order size and a smaller order point. This difference seems to be consistent when 
different examples are computed; the reason for it can be inferred from the cost function (6.2) 
and the nature of the "classical" approximations. In cases where quantity discounts are given 
for large orders, the difference may be even more pronounced. + 

In practice it is important that the cost rate be relatively insensitive to errors of 
estimate of parameters. Suppose that the policy-maker has available only estimates of the 
various cost and demand parameters and must establish R and r using them. For a very 
simple example, suppose that the policy-maker estimates the parameter values adopted above 
without error, except for the value of \. In Case A, A = 1.2x 10° (true); \ = 1.0x 10° 
(estimated). In Case B, A = 0.8 x 10° (true); A = 1.0x 10° (estimated). If the estimated value 
of X were assumed to be the true value, the policies given by (6.9) and (6.11) would apply. 
The corresponding true cost rates would be: 


Case A: C(200, 183) 89,404 (dollars per year) 
C (290,140) = 84,845 (dollars per year) 
C (200,183) = 63,949 (dollars per year) 
C (290,140) = 61,647 (dollars per year). 


In both of these cases, mis-estimation of the demand rate by approximately 20 percent, above 
and below, results in higher cost rates when the "classical" formulas are used than when the 
policy (6.11) is used. Thus (6.11) is probably relatively less sensitive to demand rate errors 
than is (6.9). In general, it is likely that a policy that is more nearly optimal than the "'clas- 
sical" policy when parameters are known will also be relatively less sensitive to errors of 
estimate of parameters. Such a property could have considerable practical importance. 
Finally, a general remark is in order. The two-bin method of inventory control is one 
of the oldest in common industrial use. In spite of its long history, however, some care 
should be exercised in its use. Care is particularly needed when the two-bin system is being 
considered to control stocks of a number of items, all of which are manufactured by the same 
production facility. The reason for this is that, in the interests of "set-up" or ordering cost 
economy, relatively large but infrequent orders for replenishment stock tend to be placed with 
the production facility. Because of chance variability in usage for the various items, these 
orders appear in a random (perhaps Poissonian) fashion. The result is likely to be a "feast or 
famine" situation at the production facility—either the facility is idle or it has a large backlog. 


; tit has been pointed out by a referee that an iterative procedure can be used to improve the 
simple "classical" result (6.10). 
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In the latter case there is temptation to use costly overtime production. The larger the set-up 
costs are made, as compared with carrying costs, the larger become the "optimal" values of 
R and the more infrequently orders arrive; but the more the production load tends to fluc- 
tuate—usually an uneconomical situation. The magnitude of these fluctuations can be estimated 
with the aid of waiting-line theory, but the need is for a system designed to reduce them to a 
tolerable level. 
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APPENDIX I 
Proof that the Inventory Distribution Attains Stationarity 


The simplest proof of the existence of a stationary distribution depends upon a 
Tauberian theorem for Laplace-Stieltjes transforms. We refer the reader to Widder [6, 
p. 192] for statement and proof of this theorem. 

To apply the theorem, take Laplace transforms with respect to the time variable on 
both sides of Eq. (3.8). By use of the convolution theorem for L. transforms, 


(I.1) F(i,s) = G(i,s) + G(i,s) Q(s), 


where 
tee) 
F(i, s) = l et Fi, t) dt, 


@ 
G(i, s) = j «= G(i, t) dt, 


© 


Qi, s) = { e* agit). 
Replacing Q(s) by its transform in terms of the basic distribution Q, (s), we have from (3.4) 


Q, (s) _ G(i, s) 
1-Q,(s) 1 - Q(s)° 





(1.2) F(i,s) = G(i, s) + G(i, s) 


To apply the Tauberian theorem, multiply both sides of (1.2) by s and take the limit as s>0. 
By the theorem, 


lim s F(i, s) = lim 3 { F (i, t) dt 
s—>0 t> 2» 


i 0) 


oe f G(i, t) dt 
ie ee 
s->0 {1 - Q,(s) co 
—sz J | taqe 


Notice that this actually proves only that a Césaro mean type of stationary distribution 


exists. The existence of the above limit (I.3) does not prove that lim F(i, t) exists, as many 
—s to 





examples show. However, if the latter limit can be shown to exist, the two must be equal. That 
the "time limit" 
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(1.4) lim F(i,t) = F(i) 
to 


exists follows from a basic result in renewal theory [3], which states that if Q: (x), the 
probability distribution of reorder times, is absolutely continuous, then 


_ dQ a olin 
(1.5) Ree _ (t) = q(t) = EC’ 


or, speaking figuratively, the "reorder density" tends to a constant. The hypothesis that 
Q) (x) be absolutely continuous (possess a density function) is a mathematical essential for 
this result to hold, but from the point of view of applications it is unrestrictive. 

Blackwell's result can be applied directly to equation (3.8) to produce (3.9). The steps 
in the argument are simple and will be omitted. The limit (3.9) is independent of initial con- 
ditions, for we let Q (t) be the distribution function of the time to pass from I = ip; the initial 
inventory, to the first cut of the order point. ‘Since Q (t) is a bona fide distribution function 
allotting no probability to infinite times, then with probability one the order point is cut at 
some finite time, and the previous analysis applies. ’ 


APPENDIX II 
The Laplace-Stieltjes Transform of the Passage Time Distribution; 
Derivation of the Average Passage Time 


If v is a fixed quantity of stock (v > 0) and T,, is the time taken for demands to first 
exceed v, we have from (4.6) 


Prob {Ty > t} = H(t, v), 
so 


(11.1) Prob {Ty < t} =z 1- Ho(t, v) = G(t, v). 


The Laplace-Stieltjes transform of this distribution function is 
i ¢) 


(11.2) G(s, v) J est d, G(t, v) 


@ 
(11.3) -G(0, v) +s f e"St Git, v) dt. 


The integrals converging at least for Re(s) >-.’. For v > 0, G(0, v) = 0, so we need only 
evaluate the integral 
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© @ 1 co 
(0.4) | et Git, v) dt = I et [1- Hy(t, v)] dt = -- I eS H(t, v) dt. 


We obtain, by referring to (4.5), 


ie 8) 
(11.5) if est Ho(t, v) dt 


1h 


This follows upon reversal of the order of integration and summation. By substituting into 
(1.4) we have 


(11.6) 


This is the transform we are after. 
Differentiating (II.6) generates the moments of the passage time. Thus 
ak G 


© 
G_ 7.4) k .-st 
djk ( 1) { te dy Git, v) 


exists by the analytic character of the Laplace-Stieltjes transform. Thus 


x a&G 


By, = (-1) af 


and we have 








THE WATCHDOG AND THE BURGLAR* 


S. P. Thompson and A. J. Ziffer 


Naval Research Laboratory 


There is presented a purposefully oversimplified game-theoretic analysis of a com- 
petitive situation which is possibly of considerable practical importance to certain engineering 
developments. For purposes of exposition, the work will be interpreted in terms of watchdogs 
and burglars. Mathematically, it is a consideration of the dependence upon two parameters of 
the value of a certain matrix game played between a watchdog and a burglar. 

These parameters are interpreted as measuring the biting ability and the stealth of the 
watchdog, and our principal interest is in determining the circumstances under which a breeder 
of watchdogs should lessen the one in order to increase the other, and vice versa. The game is 
introduced as a mathematical device which allows the determination of the optimal behavior in 
the maneuvering of the dog and the burglar for all values of the parameters. Assuming such 
optimal behavior, the payoff to the breeder then becomes a function only of the parameters, 
making possible quantitative statements about the relative desirability of different portions of 
the space of these parameters. , 

Hence consider the following situation. A burglar has intruded into the watchdog's yard. 

The dog positions himself squarely behind the burglar so that by making a last forward lunge he 
may bite the burglar with sufficient severity to render him unfit to pursue his occupation. At 
this juncture, a two-person zero-sum game is played. We suppose that the watchdog can lunge 
straight ahead, bear to the left, or bear to the right; and that the burglar is capable independ- 
ently of similar maneuvers. Also let us suppose that the burglar detects the watchdog with 
probability p. If he detects the watchdog, then he can choose any one of the three directions. 
If he does not detect the watchdog (which happens with probability 1 - p), then we assume that 
he goes straight ahead. In any case, both the watchdog and the burglar know that the watchdog 
is ignorant of whether or not the burglar has detected him. For convenience, we shall hence- 
forth refer to the eventuality of the burglar detecting the watchdog as detection. 

We express the payoff of the game as the probability of the burglar's receiving an 
incapacitating bite and, for simplicity, assumet that this probability is unity if both choose 
identical directions and the interception is therefore squarely made; is zero if the burglar 
dodges left/right when the dog lunges right/left and the interception fails; and has some 


intermediate value \ for grazing interceptions, i.e., when one continues straight and the other 
deviates in either direction. 


*Manuscript received June 13, 1958. 
At the end of the paper some gereralization of these assumptions is discussed. 
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We formalize these considerations by employing the numbers 1, 2, 3, to indicate the 
directions left, straight, and right, respectively, and by employing W for watchdog and B for 
burglar. The following array then summarizes the foregoing pay-off assignments: 


























It is to be noted that the elements may be interpreted as payments by B to W: if the numbers 
representing the directions chosen differ by two, then W receives zero; if the numbers are the 
same, then W receives unity; and in the case where they differ by one, W receives the inter- 
mediate reward of X, where \ is some number such that 0<. <1. The parameter \, the 
pay-off to the watchdog for a grazing interception, can be considered directly proportional to 
biting ability. Finally, we remark that both W and B are aware of the values of the parame- 
ters p and ) in any particular play of the game. 

The preceding description of the game is in extensive form. To put the game in normal 


form, we let Wis i= 1,2,3, be the probability of W employing his i-th pure strategy, defined 
as: 





W's i-th pure strategy: turn in direction i; and let b,, j = 1,2,3 be the probability of 
B emplying his j-th pure strategy, defined as: 


B's j-th pure strategy: turn in direction j if detection occurs, otherwise take 
direction 2. The normal form of the game may then be represented by: 





by by 





Wy |} Pir Pie 
Wo | Poi Poo 


W3 | P3i Ps2 P33 » 





where the p's are functions of p and of the elements of the array previously displayed, and 
where each represents the pay-off which would occur if the corresponding values of w and b 
were unity and the remaining w's and b's were zero. 

Consider p,,. If W lunges left then, with probability (1 - p), B fails to see W, goes 


straight, and the pay-off is \; while with probability p he sees W, goes left, and the pay-off 
is unity. Therefore 


Pyz = (L-p)ArA+p. 
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Similar reasoning shows that: 


The expected pay-off is then: 


i 


(4) Zz Pij Wy 2; » 


i,j=1 


which we assume W desires to maximize by choosing the w;'s subject only to the restriction 
that they be nonnegative and sum to unity, and which we assume B desires to minimize by 
choosing the b.'s subject to a similar restriction. The optimal strategies for W and B are 
mathematically those probability vectors w = (wy, Wo, W3) and b= (b,, bo, bg) which minimax 
the expected pay-off in the manner customary in game theory and which may be interpreted as 
indicating the probabilities that each in self-interest should employ. 

All of the optimal strategies have been obtained by examining, in the manner given by 
McKinsey [1], all square sub-matrices of (2). The character of the solutions depends upon the 
location of the parameter point (A, p) in the plane of these parameters. This plane naturally 
divides into the three regions sketched in Figure 1 and specified by: 


1 
I. P <5? 


1 2- 3p 
0. — 2 
p>—> > 301 = 2p)’ 


Ia. ose, i 
2 2(1 - 2p) 


The optimal strategies corresponding to each region are found to be: 

In I, w= (0,1,0), b= (r,0,1-r), where r is any nonnegative number not greater than 
unity. Left and right turns are superfluous pure strategies for W, as is the choice of the 
straight-ahead direction in the case of detection for B. Hence W should always go straight, 
and every mixed strategy which avoids going straight in the case of detection is an optimal 
Strategy for B. 
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Figure 1 
In 0, w= (0,1,0), b= (s,0,1-s), where s is any number such that 


(1 - 2p) (A - 1) | Pe Be. 
p orien p : 








The pure strategies which are superfluous in I are also superfluous in II. However, B can 
no longer play optimally merely by avoiding going straight in the case of detection. 

In Ill, w= (t, 1-2t,t), b= G 1- a) eve tee. 

p pp 3-42 
fluous for either player in the whole region, and the optimal strategies for both are symmetric 
in the sense that both must play their first and third pure strategies with equal probability. 

We omit from discussion the problems facing the instructors of watchdogs and burglars 
in teaching their pupils to randomize their jumps in accordance with optimal strategies. 
Assuming this to be done, we concentrate upon the problems facing the breeder of watchdogs, 
examining in greater detail the implications of the value of the game, defined in game theory 


as the expected payoff (4), evaluated under the optimal strategies. Using V to symbolize the 
value, one finds that: 


InI and, V=pA+(1-p); 


- No pure strategy is super- 


In I, y ob 22%) 
(3 - 4A) 


As is always the case, the value (interpreted as the probability of an incapacitating bite) is a 
unique function of the parameters, and inspection will show that the value is continuous across 
2-3 
the line A = en separating regions I and II (where V is a function both of \ and p) from 
- 4p 


Ill (where V is a function only of ). 
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Figure 2 


Figure 2 will assist in the discussion of these functions and their implications for the 
breeder of watchdogs. The heavy lines are isovalues (curves of constant V) in the (\,p) plane. 
The isovalues for V < 3 become of infinite slope as they cross the dotted curve separating 


region III from the rest of the unit square, since above this curve V is a function of ) but not 
of p. We recall that biting ability increases with \, while stealth decreases with p, and that 
V is the probability of an incapacitating bite. As might be anticipated, V increases with ) at 
constant p. However, watchdog breeders rarely get something for nothing, and it is to be 
expected that in practical cases the breeder will be bound by an additional relationship con- 
necting A and p whichwillprohibit him from keeping p constant as he increases \. Because 
increased biting ability in the watchdog is almost certainly obtained at the price of his easier 
detection by the burglar, because of the dog's larger size, it is to be expected that in practice 
p is a monotone nondecreasing function of ’. The given plot of isovalues may be used to 
explore the consequences of any such relationship which a practical watchdog breeder might 
desire to hypothesize. ; 

Some general tendencies controlled by the form of such relationships may be easily 
illustrated. If the function is convex downward, like curve 1 in Figure 2, it tends to be parallel 
to the isovalues themselves, and the value of the game (the probability of an incapacitating 
bite) hence tends to vary but little as one moves along the curve. In Figure 3, the values of 
the game which one obtains in this motion are plotted as a function of \, here, also, curve 1. 
While the breeder should choose ) as large as possible in this particular case, V for \ = 0.8 
is only about 35 percent larger than it is for \ = 0, and it is not obvious that the breeder should 
be greatly concerned about optimizing the breed. 
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Figure 3 


However, the variation in V tends to be larger when p as a function of ) is convex 
upward like curve 2 in Figure 2. The corresponding plot in Figure 3 now shows that V for 
\ = 0.1 is about 65 percent larger than for \ = 0.4, where it is a minimum. 

The tendency for large variation in V along the p(d) curve is even more accentuated 
for S-shaped curves like the one labelled 3 in Figure 2. The corresponding curve similarly 
labelled in Figure 3 shows that V at A = 0.1 is now about 2-1/2 times greater than its minimum 
value at 4 = 0.25. This circumstance would worry the authors if they were watchdog breeders, 
for they suspect that this is actually the type of curve to be anticipated in practice. They 
would expect as burglars to see small dogs with a low probability nearly independent of the size 
of the dog, and would expect this probability to increase very rapidly once the dog exceeded a 
critical size, saturating finally at a value less than unity, a value determined not by the dog, 
but by the burglar's imperfect alertness. (It can, of course, be true that this conjecture is 
correct, but that the saturation value of p is small enough for V to remain large and relatively 
constant over the whole of the p(X) curve, justifying less concern in selecting the value of 
than appears to be appropriate in the hypothetical example cited.) 

The deep minimum in curve 3 of Figure 3 illustrates, for the similarly labelled 
S-shaped curve in Figure 2, the danger of uncritical compromise between stealth and biting 
ability. Narrow-minded concentration on either of these attributes to the exclusion of the other 
would, in the absence of any analysis, give a watchdog of high performance, while well- 
intentioned but unthinking compromise is liable to result in unfortunate consequences. 

It is interesting to consider two generalizations of this simple watchdog-burglar game. 
The first one involves only the following change: when it comes time for W to play, we sup- 
pose that with probability © he receives a pay-off of 1, and with probability 1-o he plays as 
before. This variation is of particular interest for the case # = 1, since we may interpret the 
new feature as incorporating the probability 0 that W perceives B's maneuver and rushes the 
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interception squarely and the probability (1-c) that W does not so observe and is forced to 
play as before. 
The strategy pay-off matrix A, of the new game takes the form 


A, = (l-o)A +uo0Jd, 
where A is the strategy pay-off matrix (2) of the original game, and J is a matrix of all ones 
and of the same order as A. Owing to the fact that the elements of A, are formed by multi- 
plying all the elements of A by a common factor and then adding a common number to all, the 
optimal strategies for the A, game are the same as the optimal strategies for the A game, 
and the values of the games V(A,) and V(A) are related by the equation 
V(A,) = (1-0) V(A) + wo. 

Another generalization is a continuous analog of the original game. Let w and b be 
continuous variables in the interval 0 < w, b <1, and consider the function M formed as the 
sum of a function M, applying with probability p and a function My applying with probability 
(1-p): 


M (w, b) al p M, (w, b) + (1 -p) Mo (w, b) , 


M, (w,b) = (2 - 42) (w-b)? + (40 - 3)|w-b/ +1, 


Mg (w, b) = 4(r - 1) w? + 4(1-A)w+d. 


The values w,b = 0,5 1 are intended to correspond respectively to the maneuvers left, center, 


right in the original game. That these identifications are correct may be seen by computing 
the array 


M, (0,0) M,(0,5) Mj (0,1) 


1 11 1 
M,(5,9) M,(5,5) My (5,1) 








1 
M, (1, 0) M, (1,5) M, (1, 1) 


1 
and verifying that the result is (1); and computing Mo (0,5), Mo 9) My (1,5) and verifying 


that these constitute the center column of (1), that is, the pay-offs to W when B goes straight. 
Strategies for a continuous game consist of probability distribution functions which for 
our game we denote W(w) and B(b). We have found optimal strategies in those regions where 
there were superfluous pure strategies at every parameter point, namely in regions I and II. 
Using the step functions I, (z) and Ip (z) defined for arbitrary 0 <z <1 and 0<a<1 as 
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1, (2) when z <a, 
1, (2) when z>a, 
Ip (z) = when z=0, 
Ip (z) = when z>0, 


we find the following to be optimal strategies in regions I and II: 


W (w) - I, (w) , 
2 


B(b) = ; {Ip (b) + 1, (0)| . 


1 1 
It will be noted that these are analogs, respectively, to the vectors w = (0,1,0) and b= (5,0,5), 


which in the original game are optimal strategies in the aforementioned regions. — 

In these regions the value V(,p) in the continuous game turns out to be the same as 
the value in the original game. There are therefore grounds for belief that the previously 
discussed implications for the breeder of watchdogs stem primarily from the fundamental 
nature of the competitive situation treated and not from the fact that, for convenience, principal 
attention was devoted to a discrete formalization of this situation. 
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The Standard described here was prepared to meet a growing need 
for the use of standard sampling plans for inspecting products procured 
by the Government, It may serve asan alternate to MIL-STD-105A, which 
provides sampling plans for lot-acceptance inspectionby attributes. The 
variables sampling plans apply to a single quality characteristic and 
may be used whenever the quality characteristic for test can be measured 
on a continuous scale. Like the pattern of MIL-STD-105A, risks of 
rejecting acceptable material were planned to decrease with increase 
in size and lot. In the new Standard, however, samples sizes are con- 
siderably smallerthan those in MIL-STD-105A for corresponding sam- 
pling plans. 


The Standard outlines procedures to be followed whenthe variability 
is unknown--either a standard-deviation or anaverage-range method-as 
well as those to be used when the variability is known. A new concept 
in the Standard is a lot-acceptance criterion obtained by comparing an 
unbiased minimum -variance estimate, r, of the true percent defective with 
a suitably chosen constant, M, which represents the maximum allowable 
percent defective for sample estimates. This procedure, when used for 
the two-sided specification-limit case, provides operating characteris- 
tics, for given percent defective, which are essentially the same as 
those for the one-sided specification-limit case, regardless of the total 
percent defective in the lot, which consists of the sum of the percent 
defective above the upper and below the lower specification limits. For 
the one-sided specification-limit case,the acceptance criterion provides 
the identical probability of acceptance as obtained with the familiar 
X +ks statistic for suitably chosen k. Values of k also are included in 
the standard, thus permitting use of either form of the acceptance 
criterion. 











INTRODUCTION 


Acceptance inspection may be considered a procedure for making decisions as to 
whether to accept or reject a quantity of product submitted by a supplier. Statistical sampling 
plans have been made a vital part of such procedures, because they are an effective tool for 
providing, in advance, any desired control over the correctness of these decisions. When a lot 
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of any given quality is submitted for acceptance, and a random sample is selected, the 
operating characteristics of the sampling plan give the probability associated with the decision 
to accept the lot. 

For quality measured in terms of percent defective, inspection by attributes, where an 
item is classified as either defective or nondefective, based on inspection of specified quality 
characteristics, is most often appropriate. In this case MIL-STD-105A [1] provides an ample 
variety of sampling plans from which the inspector may select, according to definite proce- 
dures, the sampling plan required for his inspection assignment. To determine lot acceptance 
under MIL-STD-105A, the number of defective units in the sample are counted and the total is 
compared with an acceptance and rejection number. 

However, if the inspection or test of a quality characteristic can provide measurements 
on a continuous scale, and if the distribution of the measurements is known, an acceptance 
inspection procedure can be devised that will give the same control of decisions as that 
afforded by the attributes sampling plan, but at considerable reduction in sample size. Such 
measurements are often referred to as variables data, and the acceptance inspection proce- 
dure is known as inspection by variables. 


ORIGIN OF MILITARY STANDARD 414* 

In the Department of Defense, use of inspection by variables as an acceptance sampling 
procedure in order to evaluate the quality of product submitted by a supplier has increased 
markedly in recent years. During the same period, workers in the field of quality control in 
industry, at universities, and in the Government have recognized the importance of this 
sampling procedure. In the published and unpublished literature on quality control, there is 
ample evidence of the contributions to statistical research in support of this test procedure and 
of the numerous applications by various activities. Some activities in and out of the Govern- 
ment had already prepared sets of standard sampling plans designed to their own preferred 
methods of application. The difference in methods of application was undoubtedly influenced by 
the availability of theoretical results at the time of preparation or by different opinions as to 
the ease of application. Nevertheless, these differences are all superficial because, basically, 
all percent-defective inspection by variables sampling plans worthy of the name are cut from 
the same cloth. In view of this, the need for a military standard for inspection by variables 
for percent defective, to serve as an alternative to MIL-STD-105A for inspection by attributes, 
became apparent, and its preparation was a natural consequence. 

Accordingly, the Office of the Assistant Secretary of Defense (Supply and Logistics) 
brought together, as a team, representatives of the various agencies of the Department of 
Defense to specify the scope and content of the new military standard. This specification 
included types of sampling plans, lot sizes, sample sizes, acceptable quality levels (AQL's), 
risks of incorrect decision for AQL products, and the relationship among the various types of 
sampling plans. 

A desire of this task team was to use the theoretical solution to the problem of specify- 
ing an acceptance region for the two-sided test in inspection by variables with properties 
comparable with that for the one-sided test, which was first prepared by Bowker and Goode [2]. 
This led to a change in the acceptance criteria from a graphical form to a simple analytical 





*Military Standard 414, "Sampling Procedures and Tables for Inspection by Variables for 
Percent Defective,'' is available from the Superintendent of Documents, Government Printing 
Office. 
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procedure, a detailed discussion of which will appear later in the paper. The theoretical 
solution and the revision of the acceptance criteria was accomplished by members of the 
Applied Mathematics and Statistics Laboratory, Stanford University, under a quality control 
research program administered by the Statistics Branch, Office of Naval Research, and sup- 
ported jointly by the three Military Services of the Department of Defense. The Stanford group 
also completed the sampling procedures and tables within the framework indicated by the 
Department of Defense Task Team in the form of a first draft and is also responsible for the 
theoretical basis of the Standard. This work was reported by Lieberman and Resnikoff [3]. 

The Bureau of Ordnance, Department of the Navy, was responsible for preparing and 
coordinating MIL-STD-414 in a form acceptable to the Army, Navy, and Air Force. In the 
process of preparing the Standard, several revisions, incorporating certain corrections and 
additions, were made to the first draft, leading to the present format. This work was accom- 
plished by the authors while they were members of the Statistics Branch, Quality Control 
Division, Bureau of Ordnance. In addition, it was also considered desirable to publish, for 
ready reference, a technical report on the theoretical basis of the procedures and tables in 
the Standard. This report [4] has now been issued by the Department of Defense. A listing of 
all those people whose numerous helpful comments influenced the contents of the Standard in 
its present form would be long indeed. It is fair to say, however, that they represent leading 
persons in the fields of quality control and statistics from universities, industry, and the 
Government. 


BRIEF REVIEW OF INSPECTION BY VARIABLES FOR PERCENT DEFECTIVE 
AS A TEST PROCEDURE 

If the quality of a product depends on a characteristic that can be measured on a con- 
tinuous scale, such as hardness, time, diameter, or weight, it is possible to use inspection by 
variables as a test procedure. If the measure of quality is percent defective, there is specified 
an upper specification limit or a lower specification limit—sometimes both limits—on the 
measurable characteristic. Any unit of product for which the measurement of the quality 
characteristic exceeds the upper limit and/or is less than the lower limit is a defective unit. 
The lot percent defective, say p', is then the result that would be obtained if the total number 
of defective units in the lot, multiplied by 100, were divided by the lot size. The objective of 
the test procedure applied to acceptance sampling is to select a sample from the lot and decide 
to accept or reject the lot on the basis of the collective measurements of the sample units. 

To develop a rational basis for the decision procedure, some value of lot percent 
defective is specified, say Po as an acceptable quality level (AQL), such that it is preferred to 
decide to accept the lot if the true lot percent defective does not exceed this AQL value; like- 
wise, some value of the lot percent defective, say Py where Pi > Po that represents an 
undesirable measure of quality is specified for which it is preferred to decide to reject the lot 
if the true lot percent defective is equal to or greater than this value. These preferences are 
made specific by the designation of suitable risks a and 8, where a is the risk of making an 
incorrect decision to reject when the true lot percent defective is Po and f is the risk of 
making an incorrect decision to accept when the true lot percent defective is Pj: These risks 
represent the probability of an incorrect decision or, what is equivalent, the frequency in the 
long run of making an incorrect decision. An objective of the inspection by variables test 
procedure is to provide a method for making acceptance decisions that will be controlled within 
the risks specified for quality levels Po and Pi with a minimum fixed sample. The operating 
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procedures provide the size of the sample to select and an acceptance constant. By a com- 
parison of the sample results with the acceptance constant, a decision to accept or reject the 
lot is made. The procedure makes it possible to compute the operating characteristics of the 
sampling plan, that is, the probability of acceptance, L(p'), for any value of the true lot per- 
cent defective, p'. In this Standard, the computation of these probabilities is based on the 
assumption that the random variable, which represents the quality characteristic, is normally 
distributed. The scope of application can be increased for variables which are not normally 
distributed if it is possible to transform the random variable to one which does follow the 
normal law. 

Given the two pairs of values Po» @ and Py, 8, which represent the desired properties 
of an acceptance procedure, then the appropriate sample size and acceptance constant for a 
sampling plan guaranteeing these properties can be derived by the procedure described by 
Resnikoff and Lieberman [5] or in the technical report [4] issued by the Department of Defense. 
Alternately, a sampling plan can be selected from the collection of plans in the Standard by 
following the procedures of that document. The value of Pi is not specifically mentioned in 
the Standard but is tacitly implied by the point on the operating characteristic curve corre- 
sponding to a specified value of 8 = L(p}). 


COVERAGE OF STANDARD 
The Standard is divided into four separate sections, as follows: 
Section A, General Description of Sampling Plans 
Section B, Variability Unknown—Standard Deviation Method 
Section C, Variability Unknown—Range Method 
Section D, Variability Known 

Section A, the general description of the sampling plans, will always have to be used in 
conjunction with the other sections, because it provides general concepts and definitions needed 
for sampling inspection by variables. It includes such topics as Scope of Standard, Acceptable 
Quality Levels, Submittal of Product, Lot Acceptance, Sample Selection and Operating Char- 
acteristic Curves. 

Sections B and C describe the specific procedures for application of the sampling plans 
when variability is unknown. The term variability is a general term that refers to the measure 
of dispersion of the quality characteristic for the lot or process. For most applications the 
variability is unknown and is estimated either by the sample standard deviation under Section 
B or by the average range under Section C. The range method requires, for most of the 
comparable plans, more observations than the standard deviation method, but it has the 
advantage of being simpler in the computations required. 

Section D is available if the variability is known. Substantial savings in sample size 
over unknown variability plans, for the same operating characteristics, can be made in this 
case; however, the user is reminded that the requirement of known variability is a stringent 
one. From a practical viewpoint, variability may be considered "known" if it is observed to 
remain stable over a period of time, that is, to be in "statistical control." 

Like those of MIL-STD-105A, the sampling plans are indexed by lot size, inspection 
level, and AQL. There are 17 lot-size intervals, each identified by a sample-size code letter, 
grcuped into 5 inspection levels. Sampling plans are provided for 14 AQL values, ranging 
from 0.04 to 15.0 percent. Furthermore, the risk involved in rejecting material for a given 
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AQL decreases as the lot size increases. This again follows a pattern similar to that of 
MIL-STD-105A. 

The corresponding plans in Sections B, C, and D were matched as well as possible, 
under a system of fixed sample sizes, with respect to their operating characteristic curves 
and are extremely close in most cases. Hence, for a given lot size, inspection level, and AQL, 
the user may select a plan from either Section B, C, or D and be assured of essentially the same 
probability of rejecting or accepting material for any given quality. It is important to realize the 
provision that the variability has to be known for the Section D plans. 

Each of the Sections B, C, and D consists of three parts: 

(1) Acceptance sampling plans for single specification-limit case. 

(2) Acceptance sampling plans for the double specification-limit case. 

(3) Procedures for estimation of process average and criteria for tightened and 

reduced inspection. 

The single specification-limit case is encountered when the quality requirement is 
given as either an upper or a lower limit for the characteristic measured and a corresponding 
single AQL is specified. In this case, the acceptance criterion is given in two forms. For the 
same sample size, one form provides an acceptance constant k similar to that used in the 
X+ks procedure; the other form provides a maximum allowable percent defective M and 
requires estimates of the lot percent defective. Either form may be used, because for the 
single specification-limit case they provide the identical operating characteristics. 

Under the double specification-limit case, the quality requirement is given by an upper 
and a lower specification limit for the characteristic measured. If a separate AQL is assigned 
to each limit with no specific control over the total percent defective outside both limits, then 
the procedures for a single limit case can be applied to each limit separately. This is not 
ordinarily considered to be a double specification-limit case. 

In most applications of a double specification limit, one AQL is assigned to the combined 
percent defective above the upper limit and below the lower limit; in certain cases separate 
AQL values are assigned to the upper and the lower limit and some control is specified over 
the total percent defective. Procedures for both the single and the separate AQL situation for 
the double specification-limit case are provided in each of the Sections B, C, and D. Additional 
discussion of the problems encountered under this case is given later in the paper. 

Estimation of process average under procedures of the Standard is accomplished by 
simply averaging the estimated lot percent defectives for a number of lots. Procedures for 
tightened and reduced inspection, based on inspection results of the previous 5, 10, or 15 lots, 
also are included. 


DISTINGUISHING FEATURES 

A distinguishing feature of the Standard is the introduction of acceptance criteria 
obtained by comparing the estimate of the percent defective p, as computed from the sample 
results, with a maximum allowable percent defective M obtained from the tables in the Stand- 
ard. Lot disposition then is as follows: 

If p< M, accept the lot. 

p> M, reject the lot. 

The maximum allowable percent defective M is an acceptance constant that depends on the lot 
size, inspection level, and applicable AQL. In the single specification-limit case, p is an 
estimate of the percent defective above the upper or below the lower specification limits. 
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This acceptance procedure is analogous to attribute inspection under MIL-STD-105A, 
where the number of defectives found in the sample is compared with the acceptance number, 
A comparison is set up in Chart I. 


CHART I 
Comparison of Attribute Inspection with Variables Inspection 
(Lot size, 3000; inspection level, I; AQL, 1 percent) 


Attribute Inspection Variables Inspection 
Under MIL-STD-105A Under MIL-STD-414 





Sample size 
code letter L 


Sampling plan Variability Unknown 
Standard Deviation 
Method 


Sample size n= 150 n= 40 


Acceptance 


a 
oti ant ( = 2.67%) M = 2.71% 


Sample results used d = number of defec- p=estimate of lot percent 
tives in sample defective based on X 
and s 


Acceptance criteria Compare d and a Compare p and M 


d 
Lot disposition: Accept Ifd<4 (< <2.61%) If p < 2.71% 


d 
Reject Ifd>4 (< > 2.67% | If p > 2.71% 


The percent defective in the lot is estimated by computing a quality index, Q, from the sample 
results. As an example, the quality index for the variability unknown-standard deviation 
method and a single specification limit, U, is: 
uU-xX 

Q= — 
where X and s are the sample mean and standard deviation, respectively. The value Q is 
then used to read off the estimated percent defective, p, from a table provided in the Standard. 
The estimate p obtained in this manner is the best possible in the sense that it is the unique, 
unbiased, minimum-variance estimate of the lot percent defective. The estimate p obtained 


by using the known standard deviation method also is an unbiased, minimum-variance estimate 
of the lot percent defective. 
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100 


PERCENT OF LOTS EXPECTED TO 
BE ACCEPTED 


| 


MEAN 











a 
rae ss 
6 


Figure 1 - Example of double speci- 
fication limit (normal distribution of QUALITY OF SUBMITTED LOTS IN 
items) PERCENT DEFECTIVE (p') 


U = Upper specification limit Figure 2 - Typical operating-charac- 


teristic curve (AQL, 1%; lot size, 3000; 


Lower specification limit sample-size code letter, L) 


Percent of items above U Variables sampling plan, 


standard-deviation method 
(Sample size (n) = 40; maxi- 
‘ ; } mum allowable percent 
= Py + » * percent of items out- defective (M) = 2.71%) 

side of limits 


Percent of items below L 


Attribute sampling plan, 
single sampling 

: : (Sample size (n) = 150; maxi- 
Most important, however, in sam- “am stein eumthae of 
pling inspection by variables is the fact that defectives (a) = 4) 
the use of this estimate p provides a solution 
to the double specification-limit problem 
when a single AQL is specified. This problem consists in providing an acceptance procedure 
for the double specification-limit case so that the operating characteristics are the same as 
those for the single specification-limit case. As an illustration of what is involved, Figure 1 
shows a normal distribution curve with the upper limit indicated by U and the lower limit by L. 
The lot percent defective above U and below L are indicated by Py and Pi, respectively. The 
mathematical procedures necessary to compute the exact operating characteristic (OC) curve 
require knowledge of Py and Pr, separately and not their sum. Hence, there is a separate OC 
curve for each split-up of the total percent defective p' which equals Py + Py: The procedures 
of the Standard, through the use of the unbiased minimum variance estimate p, are such that, 
for practical purposes, almost identical OC curves are obtained for a fixed value of p', regard- 
less of this split. No other known procedure has this property with such exactitude. The oper- 
ating characteristic (OC) curve of a sampling plan is shown in Figure 2. Similar curves are 
given for all the plans in the Standard. It should be noted that the curves are based on the 
assumption that the measured characteristic of the individual items follow a normal distribu- 
tion. The curves represent the percent of lots expected to be accepted for a given percent 
defective. Chart II gives a detailed example of the procedure taken from the Standard for a 


double specification-limit case in which the variability unknown-standard deviation method is 
used, 
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CHART II 
Example of Calculations for Double Specification-Limit Case 
by the Variability Unknown-Standard Deviation Method 
(One AQL value for both upper and lower specification limit combined) 


Example: The minimum of 180°F and the maximum of 209°F is specified as the 
operating temperature for a certain device. A lot of 40 items is submitted for 
inspection. Inspection level I, normal inspection with AQL = 1% is to be used. 
From the Standard, it is seen that a sample of size 5 is required. The measure- 
ments obtained are 197°, 188°, 184°, 205°, and 201°. Disposition of the lot is to 
be determined. 


Information Needed Values Obtained Explanation 








Sample Size: n 5 
Sum of Measurements: » X 975 


Sum of Squared Measurements: 1 x? 190,435 


Correction Factor (CF): (2X)2/n 190,125 (975)2/5 


Corrected Sum of Squares (SS): 1 x2cF 310 190,435-190,125 
Variance (V): SS/(n-1) 77.5 310/4 

Estimate of Lot Standard Deviation(s): VV 8.80 V77.5 

Sample Mean (X): © X/n 195 975/5 

Upper Specification Limit: U 209 

Lower Specification Limit: L 180 

Quality Index: Qy = (U-X)/s 1.59 (209-195)/8.80 
Quality Index: Q. = (X-L)/s 1.70 (195-180)/8.80 
Est. of Lot Percent Def. above U (Py) 2.19% From Standard 
Est. of Lot Percent Def. below L (P;) 66% From Standard 


Total Est. Percent Def. in Lot (p): 
P= Py + PL 2.85% 2.19% + .66% 


Max. Allowable Percent Def. (M) 3.32% From Standard 


Acceptance Criterion: Compare 
P= Py + Py, with M 2.85% (3.32%) From Standard 


Lot Disposition: Accept the lot, because p = Py + Py, is less than M. 
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An alternate procedure for the two-sided test had been suggested by the authors before 
the theoretical basis for the present two-sided test was developed. This procedure provided a 
maximum allowable standard deviation (MASD) in conjunction with two one-sided tests as the 
acceptance criteria and resulted in an adequate approximation [6] to the one-sided test for 
many applications. The procedure was incorporated in Naval Ordnance Standard 80 [7]. It was 
also made a part of the acceptance criteria in a sampling standard issued by the Army Ord- 
nance Corps [8]. The MASD has now served its original purpose, however, and is no longer 
necessary as part of an acceptance criteria, in view of the superiority of the two-sided test 
procedure incorporated in MIL-STD-414. Both the Navy Bureau of Ordnance and Army Ord- 
nance Corps standards are now replaced by the new Military Standard. Nevertheless, the idea 
of a maximum standard deviation (MSD) was considered by some to have value to users of the 
new Standard in indicating an upper bound for sample values of the standard deviation, and 
therefore a procedure for its computation has been incorporated in the Standard. It is thought 
that sample values of the standard deviation which are less than the MSD help to ensure, but 
do not guarantee, lot acceptance. 

For the future one can envision extension of the procedures of the new Standard for 
combining the results of several characteristics inspected by variables and also for combining 
percent defectives obtained from attribute and variables inspection. Research for this purpose 
is presently under way in the quality control research program previously mentioned. 


INTERPRETATION OF LOT ACTION 
Decisions as to acceptance or rejection made as a result of applying the acceptance 


criteria of variables sampling plans pertain to the magnitude of the percent, p', of measure- 
ments from some unknown normal population that fall outside specification limits. These 
decisions affect the lot from which the sample was drawn. The fact that the size of the lot is 
actually finite, whereas the theoretical range of a normally distributed random variable is 
infinite, sometimes raises a question as to interpretation of lot action. The interpretation to 
be placed on lot action for the sampling plans in MIL-STD-414, and the considerations on which 
it is based, is the following. 


The manufacturing process, from which the lot was formed, potentially or conceptually 
at least, can produce an infinite quantity of product. Strictly speaking, when the population is 
this infinite quantity, then the acceptance decisions pertain to the quality of the manufacturing 
process potential. If we assume that the average percent defective over many lots is essentially 
equivalent to the p' of the process, we can say that lot decisions are made with respect to 
average lot quality. Practically speaking, this interpretation of lot action should be kept in 
mind in those cases where the lot size is small, especially where there is a continuous supply 
of product in lots. As such, it is not an unreasonable basis for lot action, since an objective of 
lot rejection is to encourage the manufacturer to improve his process. Where the lot size is 
large it can be considered that the lot itself is equivalent to the population and that acceptance 
decisions pertain to the percent defective in the lot. This last consideration also requires that 
the sample size be small relative to the lot size, since their relative magnitude is important in 
another closely related question. This question is concerned with the effect of finite lot size on 
the probability of lot acceptance. 

The practical magnitude of "small" and "large" for lot sizes, in this context, has not 
been fully developed in the literature for the case of variables inspection. No definitive rule 
may be given without further clarification from exact studies which are needed to answer these 
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questions. It can be expected, however, that where “normality” is obtained within a finite lot, 
(in the sense of a normal bowl of chips) and sample units are selected at random without 
replacement, the risks, say a* and 8* actually obtained for the finite population (the lot) are 
never greater than the corresponding risks qa and £, respectively, for an infinite population, 
and that the difference is negligible where the sample size is small relative to the lot size. 
This implies that where normality is obtained within the lot itself, the operating characteristic 
curve with respect to inferences on lot quality is more discriminating than the operating char- 
acteristic curve given in the Standard with respect to inferences on process quality. Obviously, 
if for a suitably chosen random variable, the normality assumption cannot be made for the 
infinite population or for finite lots, particularly sufficiently small lots, any OC curve for the 
infinite population or for the finite lot based on an assumption of normality is incorrect, and 
the correct OC curve can only be computed, by using the applicable distribution, if it is known, 
in each case. In attribute inspection, the exact distribution for finite lots is known, and hence 
the exact probability of lot acceptance can be computed. When this is done decisions as to 
acceptance or rejection pertain to the quality of the lot itself and not that of a larger popula- 
tion. However, in inspection by attributes using MIL-STD-105A 1, the interpretation of lot 
action is the same as that given above for MIL-STD-414. 
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PROBLEMS 

It has been suggested that this journal might serve as a medium for publishing problems 
in the area of Logistics — with the idea that other persons might be interested in those prob- 
lems and might submit comments. Readers are invited to submit brief statements on applied 
and theoretical problems in Logistics. Address letters to Managing Editor, Naval Research 
Logistics Quarterly, Office of Naval Research, Washington 25, D. C. 





The following letter has been received: 


Dear Sir: 

The questions raised by McShane and Solomon in their letter to the N.R.L.Q. (December 
1958) concerning the Naval electronics allocation problem seem much to the point — as does 
their praise of the original papers. We should like to comment further on the potential ambi- 
guity of priority lists, such as the MIP when used to rank allocation plans. As McShane and 
Solomon have pointed out, it is not at all clear precisely what implication the priority informa- 
tion has or should have in certain concrete cases. The variety of possible interpretations of the 
meaning of such priority lists is even greater than that suggested in their letter — although they 
do not all seem equally acceptable. 

To explain, consider the following situation. Suppose there are two types of activities, 
where type 1 has been given a higher priority than type 2. Suppose further that there are three 


instances of each of the two types of activity. For present purposes assume that there is only 


one model type available. Let a hat be x 


ij instances of the assignment of a type i model toa 


type j activity in some assignment plan. Then %14214 + X97 10 specifies the number of models 
assigned to the six places in some particular plan. 

Now McShane and Solomon present two interpretations of the intent of the priority ratings 
as a source of decision-making rules for ranking assignment plans. First: “Is it the intent of 
the MIP that if it can be done, all priority 1 requirements be fulfilled before priority 2. . .?” 

(p. 363). We find it possible to interpret this first “interpretation” in at least two ways. The 
first would be: “When there is a fixed number of models available, all the higher priority 
requirements should be filled before the lower ones are even partially filled”. Thus in the 
following pair of alternative assignment plans, 


a- (3)asy + (0)aj5 


b - (2)ayy + (1)ajo 
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the rule would pick plan a over plan b. A second reading would be: “A plan is to be preferred 
to another plan if it contains a higher number of assignments of equipment to the higher priority 
activity, regardless of how many assignments to the lower priority activity are involved”. Thus, 
with a pair of plans a (above) and c, 

c- (2)asy + (2)ayo 
the rule would determine that plan a should be given a higher “military worth” rating than plan 
c. The first interpretation would be silent in ranking the two plans. 

McShane and Solomon’s second suggested interpretation of the intent of the MIP is 
phrased as follows: “Is it the intent to maximize the total ‘effectiveness’ of all the ships covered 
by the particular MIP with due consideration given to the specified ‘priorities’” ?(p. 364). Again 
suppose activity 1 has the higher priority rating and consider the following pair of possible 
assignment plans: 

d- (2)a,, + (O)ayo 


e- (1)as, + (1)ajo- 
Under the second interpretation given by the writers e could be preferred to d, while this is 
impossible under either reading of their first interpretation. 

Thus several different and conflicting rules of thumb can be deduced concerning the use 
of the MIP. Is there a basis for deciding among them [If the MIP (or any other priority list) 
should really be taken as a directive establishing a method for ranking allocation plans, then 
the natural resolution of the problem would be merely to ask the relevant officials which rule 
they wished to have applied. It would seem more reasonable, however, to regard priority lists 
simply as an ordering of importance of different activity classes — i.e., as an expression of 
military judgment regarding the ranking, not directly of alternative assignment plans, but of the 
individual ai which compose them. In this case we are left with the task of deciding which rule 
of thumb for ranking assignment plans gives results which are consistent with the expression of 
the order of importance of the activities. 

It appears to the writer that a major stumbling block in the determination of a rule for 
ranking allocation plans on the basis of priority information is the possibility of a military 
equivalent of the economist’s “diminishing marginal utility”. Indeed this is the source of the 
dilemma posed here by McShane and Solomon. “Consider the pair of assignment plans d and e 
again. Since the priority list tells us that any is regarded as of more military importance than 
aio plan e could be preferred only if we have diminishing marginal utility of pieces of equip- 
ment assigned to a given activity class. Whether this is actually the case depends upon the con- 
crete problem. If, for example, we interpret “(2)a,,” as “the installation of two radar sets on 
a single ship,” then diminishing marginal utility is probably an unavoidable implication. On the 
other hand, if the notation means “the installation of two radar sets on two distinct submarines 
of the same priority class”, then I think we might hope to rule out diminishing marginal utility 
and likewise the possibility that plan e could ever be preferred to plan d. I think it will be found 
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that j that the information contained in such priority lists as the MIP will be of comparatively 
little value in deciding among complex alternatives if we cannot assume some kind of constant 
marginal utility. 

I think that the adoption of this constant marginal utility assumption logically leads to 


the following rule: “One assignment plan shall be considered preferable to another if each 


particular aii in one has corresponding to it an aii in the other (different in each case) to which 


it is preferred or indifferent (not all indifferent). A preferred assignment plan must contain at 
least aS many ai as the plan to which it is preferred”. This rule corresponds in spirit to the 
Pareto optimality rule of economic welfare theory. 

The rule appears to have the same consequences as that found in our first reading of 
the McShane and Solomon interpretation number one — in all cases covered by the latter. It 
would be unable to decide between plans a and c, however, and it would pick d over e. As the 
case of a and c shows, it does not pick a preferred choice for all conceivable pairs of alterna- 
tives. It is, however, evidently the strongest rule obtainable so that choice between over-all 
assignment plans reflects the priority lists (interpretated as orderings of activities) and con- 
tains no additional arbitrary judgments. It should be remembered, of course, that the whole 
problem is considerably more complicated where a variety of potential model types exists. 

The writer is currently engaged in what is hoped will be a thorough investigation into 
the logical underpinnings of the use of priority lists to obtain maximizable objective functions 
in military allocation problems. We are particularly interested in the realistic case where 
only orderings are meaningful. Any communications from others who are working on or have 
thought about this problem would be greatly appreciated. 


Richard G. Davis 
Econometric Research Program 
Princeton University 





ERRATA 


List of Corrections for the paper 
“The Allocation of Development Funds: An Analytic Approach” 
by Scott C. Daubin 
Volume 5, Number 3, September 1958 


Old Version Correct Version 
Page 268, Line 16 Rate of Minimum Cost - R. Rate of Minimum Cost - oj 
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Page 268, Line 23 


Page 269, Eq. (11) 





Page 269, Eq. (12) 





Page 269, Eq. (13) 


Page 271, Line 5 


Page 271, Eq. (14) 





Page 272, Line 7 


from bottom — 





2r 


3 


. + 100r 


3 


2, 20,000 








+- 


— +— + 


2 


< 





2 
2re + 100r 


3 


+ 20,000 


+— 





NEWS AND MEMORANDA 


Readers are invited to submit to the Managing Editor items of general interest in the 
field of logistics 


For information of the members of The Institute of Management Sciences, the Sixth Annual 


International Meeting will be held in Paris in September 1959. Inquiries should be addressed to: 


Professor Wallace W. Gardner 
School of Business Administration 
University of Michigan 

Ann Arbor, Michigan. 








RECENT PUBLICATIONS 


STATISTICS FOR MANAGEMENT, by B. J. Mandel. Dangary Publishing Company, Balti- 
more, 1956. xii+ 408 pp. 

With the growing recognition of the power of statistical analysis in industrial management, 
texts for beginners, particularly those of a nontechnical nature, are in demand. To comply with 
this demand, Mr. B. J. Mandel has prepared Statistics for Management, subtitled “A Simplified 
Introduction to Statistics.” The author has achieved the aim of the subtitle, for the book con- 
tains a minimum of mathematics; the various topics are presented in simplified technical 
language. The way in which each formula is used is spelled out in complete detail. While this 
oversimplified and nonmathematical approach offers no stimulus to the more sophisticated 
student, it does offer a systematic presentation which is eash to comprehend at the elementary 
level. However, it must be pointed out that one usually pays the price for simplification. There 
are some errors and inaccuracies, particularly in the sections dealing with sampling and esti- 


mation, and some basic concepts (e.g., fiducial limits, parameters as the statistic to be estima- 


ted) are not treated at all. Furthermore, while probability is mentioned, little attempt has been 


made to explain probability as the basis for statistical theory. This reviewer believes that a 
firmer basis for statistical theory and techniques could have been established without disturbing 
the essential simplicity of the text. 

While Mr. Mandel has generally not succeeded in developing the whys and wherefores, he 
has used his illustrative material (from the field of business) remarkably well. The book is 
enhanced by the inclusion at the end of each chapter of an ample supply of pertinent problems. 
The following list of the chapter headings will give the reader an idea of the scope of the book: 
(1) The Field of Statistics; (2) Scientific Methods of Collecting Statistics; (3) Sources of Statis- 
tics; (4) Organizing Data for Analysis; (5) Averages and their Uses; (6) Variability - Its Meas- 
urement and Meaning; (7) Scientific Sampling - A shortcut to Fact- Finding; (8) Some Basic 
Principles of Sampling: (9) Estimating from Samples — Sampling Error and Risk; (10) An Intro- 
duction of Estimation by Correlation Analysis; (11) Time Series and Forecasting; (12) Index 
Numbers and their Uses; (13) Presenting Statistics in Graphs and Tables; and (14) A Summary 
of the Book. 

Statistics for Management is certainly a worthwhile text for students who wish a nontech- 


nical dose of business statistics. 


(Reviewed by P. Desind) 





STUDIES IN LINEAR AND NON-LINEAR PROGRAMMING, by K. J. Arrow, L. Hurwicz, 
and H. Uzawa. Stanford University Press, Stanford, California, 1958. 229 pp., $7.50. 

This a companion book to “Studies in the Mathematical Theory of Inventory and Production’ 
by Arrow, Karlin, and Scarf (see revies NRLQ, Vol. 5, No. 4). 

It is by no means an elementary book; those entering the field for the first time would do 
well to start with a text book such as “Linear Programming, Methods and Applications” by 
Saul I. Gass. (See review in this issue.) 

The format of the book is a series of integrated papers presenting results that have been, 
for the most part, obtained in the last three years. These papers are published here for the 
first time. 

The book deals with “fundamental existence theorems, including their extension to infinite- 
dimensional spaces as well as extension and simplification of earlier theorems for finite- 
dimensional spaces; the gradient method as a method of successive approximations for solving 
constrained maximization problems; and the exploitation of the special features of specific 
programming problems to simplify solutions.” 

There can be no doubt that the results presented will be of great value to workers in the 


field of linear programming. 


(Reviewed by J. Campbell) 


LINEAR PROGRAMMING: METHODS AND APPLICATIONS, by Saul I. Gass. McGraw- 
Hill Book Company, Inc., New York, 1958. xii+ 223 pp. 

Since 1947, when the general problem of linear programming was formulated by George B. 
Dantzig and a systematic procedure devised by him for solving it, interest in the subject has 
developed at a remarkable rate. In recognition of this interest, courses on linear programming 
are offered now in various universities, creating a need for a textbook at the advanced under- 
graduate and elementary graduate level. Although several books on linear programming have 
previously been published, none is adequate as a text. Gass’ book should fill the need in a most 
satisfactory way. It is lucid, well organized, and pleasantly written. It starts with an easily 


understood description of the mathematical nature of a linear programming problem, illustrated 


by several didactic examples (Chapter I), followed by a chapter containing those topics in 


linear algebra and convex analysis basic to the theoretical aspects of linear programming. This 
completes Part I of the book. Part II contains nine chapters on computational procedures, 
duality theory, and parametric programming. The computational methods discussed include the 
simplex method in its original, revised, and dual forms, together with a discussion of degeneracy 
procedures and a brief description of available digital-computer codes. Part III, the final sec- 
tion, deals with applications to a wide variety of problems concerning transportation, production 
scheduling, inventory control, etc., and concludes with a brief description of the connection 
between linear programming and the theory of games. Each chapter is supplemented by a short 
set of exercises. In addition, there is a valuable nine-page bibliography on linear-programming 
applications followed by a detailed list of general references. 
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The didactic value of the book is due to careful writing and skillful interlacing of theoretical, 
computational, and applied material. Nevertheless, there are occasional misleading, or at least 
ambiguous, assertions. For example, on page 56 the author perpetuates the myth that “approxi- 
mately m changes of basis are required to go from the first feasible solution to the minimum 
solution.” If this is meant as a mathematical theorem, it is false. On the other hand, if it is 
intended as a Summary of experience, that fact should have been made clear to the reader. 
Careless statements of this kind, notwithstanding, the book is generally authoritative. 

Any judgement about the selection of material is necessarily subjective. On the whole, I 
consider the author’s selection good; but, I should like to make a plea for the inclusion in future 
editions of the recent fascinating applications of linear programming to various problems in 
combinatorial analysis, such as systems of representatives and flows in networks. 

Viewing the book as a whole, we should be indebted to Mr. Gass for providing a truly 


excellent textbook in an area where none previously existed. 


(Reviewed by W. M. Hirsch) 
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